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Abstract 

Heavy-quark effective theory (HQET) is applied to lattice QCD with Wil- 
son fermions at fixed lattice spacing a. This description is possible because 
heavy-quark symmetries are respected. It is desirable because the ultravio- 
let cutoff 1/a in current numerical work and the heavy-quark mass thq are 
comparable. Effects of both short distances a and I/tuq are captured fully 
into coefficient functions, which multiply the operators of the usual HQET. 
Standard tools of HQET are used to develop heavy-quark expansions of lat- 
tice observables and, thus, to propagate heavy-quark discretization errors. 
Three explicit examples are given: namely, the mass, decay constant, and 
semileptonic form factors of heavy-light mesons. 

PACS numbers: 12.38.Gc, 12.39.Hg, 13.20.-v 
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I. INTRODUCTION 



One of the most vital parts of high-energy physics is the study of heavy quarks. Sev- 
eral large experimental data sets of hadrons with 6-flavored or charmed quarks are avail- 
able now, or will be soon. These data are valuable, because the decay properties of these 
hadrons depend on poorly known elements of the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix. A broad range of measurements can be used to determine the CKM matrix with many 
cross checks and, thus, to test the flavor structure of the standard model, including the 
origin of CP violation. 

In this enterprise numerical lattice QCD plays the role of providing hadronic matrix 
elements, ideally with controllable, transparent uncertainties. The two sources of uncertainty 
that attract the greatest concern are discretization effects and the quenched approximation, 
in which the feedback of (light) quark loops on the gluons is omitted. Both can be eliminated 
with ever larger computing resources. Effects of the lattice discretization also can be studied 
theoretically. They appear at short distances, so they can be disentangled from long-distance 
physics with field theoretic methods. This paper uses effective field theory to separate the 
short-distance scales of the lattice spacing and the heavy quark mass from the long-distance 
QCD scale. It is an extension of work with El-Khadra and Mackenzie on massive fermions 
in lattice gauge theory. The effective field theory approach yields several concrete results, 
which illustrate strategies for reducing cutoff effects of heavy quarks. Numerical calculations 
may then focus computer resources on incorporating the light quark loops. 

The idea of using effective field theory to study cutoff effects goes back to Symanzik . 
For any lattice field theory his idea was to introduce a local effective Lagrangian, which is the 
Lagrangian of the corresponding continuum field theory, augmented with higher-dimension 
operators. Coefficients of the operators depend on the underlying lattice action, and their 
dimensions are balanced by powers of the lattice spacing a, which is the only short distance 
in Symanzik's analysis. For small enough a one should be able to treat the higher-dimension 
terms as perturbations and express observables of the lattice theory as an expansion in terms 
of continuum observables. For a pedagogical introduction, see Ref . || . 

The aim of this paper is to understand numerical data generated using lattice actions 
with Wilson fermions [|]|| for the heavy quarks.Q Many papers, starting with the work of 
Gavela et al. M and Bernard et al. J7[, have attempted to calculate properties of heavy- 
quark systems in this way. In practice, the bottom quark's mass in lattice units is large, 
m^a ~ 1, and even the charmed quark's mass is not especially small, m c a « |. Thus, these 
calculations are hardly in the asymptotic regime mga — > (tjiq fixed) for which these actions 
were originally devised. In particular, any expansion in small mQa, as is usually assumed 
in analyses based on Symanzik's work, fails. This does not imply that heavy-quark cutoff 
effects in these calculations are large, but it does mean that a different analysis is needed. 

The heavy quark masses are larger than Aqcd, so they introduce additional short- distance 



In this paper the term "Wilson fermions" encompasses any action with Wilson's solution of the 
doubling problem ||. These include the Sheikholeslami-Wohlert ("clover") action ||, the actions 
of Ref. ]]J, and — of course — the Wilson action. 
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scales. One is free to seek an effective theory that lumps the effects of all short distances — the 
lattice spacing and the heavy quarks' Compton wavelengths — into coefficients. This effective 
theory does not have to be continuum QCD. The crucial observation [|1|] is that lattice actions 
with Wilson fermions satisfy the same heavy-quark symmetries || as continuum QCD. For 
heavy-light systems, therefore, a version of the heavy-quark effective theory (HQET) is 
appropriate. Similarly, for quarkonia the same Lagrangian applies, but with the power- 
counting of non-relativistic QCD (NRQCD). The operators in such a description of lattice 
gauge theory are the same as in the usual NRQCD P-pH or HQET fT2|-|15| descriptions 
of continuum QCD, but the coefficients differ because the lattice modifies the dynamics at 
short distances. 

This paper focuses on hadrons with one heavy quark and, consequently, on HQET. It 
uses tools of the usual HQET to derive formulae of the form 

Aat = 2i(m Q a) J B 00 (A QCD a) H -5^(A QCD a) H , (1.1) 

m 2 (m Q a) 

where B\ &t is a physical observable calculated in lattice gauge theory. The quantities zi(rriQa) 
and I/7712 (ttiqci) are short- distance coefficients of mass dimension and —1, respectively.^ 
They do not depend on the light degrees of freedom. The quantities B^ and B'^ describe 
the long-distance physics. They are matrix elements in the infinite-mass limit and do not 
depend on niq. Thus, the heavy-quark mass is entirely isolated into the coefficients. 

The logic to derive formulae like Eq. ( |1 . 1| ) parallels that of the standard HQET. In both 
cases the deviations from the infinite-mass limit are expressed as a series of small corrections. 
Each term consists of a short- distance coefficient multiplying a long-distance matrix element 
of the infinite- mass limit. From this structure a simple picture of cutoff effects emerges. The 
heavy-quark cutoff effects lie in the difference between the short- distance coefficient functions 
and their values in continuum QCD. On the other hand, matrix elements of the infinite- 
mass limit, such as B^ and B'^, suffer from discretization effects only of the light degrees 
of freedom. 

This paper is organized as follows: Sec. |H| clarifies the non-relativistic interpretation of 
Wilson fermions introduced in Ref. [JXJ] , by giving more direct, though also more abstract, 
reasoning to relate lattice gauge theory to HQET. Section |T| establishes some general no- 
tation and introduces the HQET Lagrangian. As in Sec. [TT] the emphasis is on symmetries. 
The leading, heavy-quark symmetric, effective Lagrangian is shown to be the same for lat- 
tice gauge theory as for continuum QCD. This static Lagrangian is the foundation of the 
heavy-quark expansion, so some of its properties are recalled in Sec. [TV]. The next task is 
to propagate deviations from the static limit to observables, so Sec. |V| develops a suitable 
form of perturbation theory. Applications of the formalism are in Sees. [VTh]VllJ| . Section [VT| 
works out the heavy-quark expansion for hadron masses to second order. Semileptonic form 
factors, at the so-called zero- recoil point, are addressed in Sec. |V11|. As with continuum 



2 Here tuq the heavy quark mass in some scheme, for example the bare mass. When there is more 
than one heavy quark in the problem, Eq. ( |f ,1| ) is schematic, and the coefficients depend on all 
heavy quark masses. 
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QCD, the first order vanishes. The technical details of the second order are considerable 
and appear in Appendix |A], correcting some minor errors in the literature. Section [VT l I 
derives the first-order expansion for decay constants. In all three cases the analysis fol- 
lows work on the usual HQET, but keeping careful track of the HQET coefficients. Some 
implications of these concrete results are discussed in Sec. IX. 



II. HQET FOR LATTICE QCD 

In this section lattice gauge theory, with a certain class of actions for Wilson fermions, is 
related to HQET. A derivation starting from the path integral of lattice QCD and making 
field redefinitions has been given in Ref. 0. That procedure is analogous to derivations of 
HQET from the path integral of continuum QCD [P?f-P0[|, and it yields the coefficients at 
the tree level. Reference [Q used heavy-quark symmetry only to show that the approach to 
the infinite-mass limit is smooth and stable in the presence of radiative corrections. Here the 
argument is reversed: owing to heavy-quark symmetry, there must be a version of HQET 
whenever momentum transfers are much smaller than itlq. Whether rriQa <C 1, rriQa 3> 1, 
or rriQa ~ 1, the reasoning is the same. The concepts are spelled out in this section, and 
the mathematical formalism is developed in Sec. |TJ. 

To show that HQET is applicable, it is enough to show that the underlying theory has a 
set of states possessing the (approximate) Isgur- Wise symmetries f| , and that no undesired 
states appear in its spectrum. Let us consider the spectrum first, because lattice field 
theories sometimes contain spurious states. The best-known are the doubler states. The 
Wilson action employs projectors |(1 ±74) to eliminate them completely |4|]: for any three- 
momentum there is only one pole in the propagator. Extra states also could arise if the action 
were to include hops over two or more time slices. Without fine tuning, the propagator then 
would have extra poles with energies near a -1 . If mqa ~ 1, as in most numerical calculations 
of heavy quarks, the effective theory would have to describe the spurious states along with 
the desired ones. This potential obstacle is easily circumvented JIJ , however, by choosing an 
action with only single hops in the time direction. 

The symmetries are revealed by writing the lattice action in the form 

S = Yl ^X^X - K ^xM x ylf)y, (2.1) 

x x,y 

where x and y run over all lattice sites and M xy has support only for y near x. To maintain 
gauge invariance M xy includes parallel transport along some path from x to y. The hopping 
parameter k controls the fermion's motion through the lattice, and small k corresponds to 
large rriQa. For k — > the propagator reduces to that of the static theory [Q, as long 
as there are only single hops in the time direction. For small, non-zero k there are spin- 
dependent and, with more than one ip field, flavor-dependent corrections. In this way, the 
approximate heavy-quark symmetries emerge naturally for large mqa. Furthermore, if M xy 
is constructed not to have unwanted poles in the quark propagator, as outlined above, on- 
shell Green functions depend smoothly on rriQa, for all rriQa > p,2~H|. Then, the Isgur- Wise 
symmetries persist as a is reduced, and one can always use a version of heavy-quark effective 
theory to describe processes with momenta small compared to rriQ. 
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Because the degrees of freedom and the symmetries are the same as in the continuum, 
the operators of this HQET are the same as those of the usual HQET describing continuum 
QCD. To define the operators the main issue is to regulate divergences. There is no need 
choose the same ultraviolet regulator for the effective theory as for the underlying theory. 
One is free to regulate the ultraviolet with, say, dimensional regularization and either a 
physical or a minimal renormalization scheme. On the other hand, because the effective 
and underlying theories are supposed to describe the same long-distance physics, the same 
infrared regulator, when needed, should be chosen. 

Because the details of the short-distance dynamics are those of the lattice theory, the 
coefficient functions of HQET must be modified. The lattice breaks some rotational and 
translational symmetries, so at non-zero a coefficients of corresponding operators need not 
vanish, as they would in the usual HQET. As a varies the short- distance properties change, 
and so the coefficients must change to compensate. Eventually, when a — > 0, lattice QCD be- 
comes (indeed, defines) continuum QCD, so the coefficients of the modified HQET smoothly 
turn into those of the usual HQET. The explicit form of the coefficients is not needed in this 
paper, but one should note that they can be calculated by computing observables in lattice 
QCD and in the modified HQET and matching. With Feynman diagrams, for example, 
one would expand lattice amplitudes around the static limit in small momentum transfers, 
keeping the full dependence on m,Qa. 

In Eq. (|2.1|) the hopping matrix M xy is not specified in detail. In general it contains many 
free couplings, which are irrelevant in the sense of the renormalization group. In the usual 
improvement program |16| they are chosen to accelerate the approach to the continuum limit. 
In the HQET analysis advocated here, a similar principle holds. The irrelevant couplings of 
the lattice action alter the short-distance coefficients of the modified HQET. Thus, they can 
be adjusted so that the HQET expansion of lattice QCD systematically reproduces more 
and more of the HQET expansion of continuum QCD. 

For a generic lattice action, the heavy-quark symmetries hold only in the rest frame. On 
a superficial glance this is a drawback, because much of the power of HQET comes from 
boosting heavy-light hadrons to arbitrary frames. On a second glance, it may be a blessing in 
disguise. By combining heavy-quark symmetry, Lorentz covariance, and reparametrization 



mvanance 



PH, it may be possible to develop a non-perturbative improvement program. 



III. NOTATION AND FORMALISM 

This section reviews the main ingredients of HQET in a notation well-suited to Euclidean 
space-time. The details are slanted to Euclidean space-time because the aim of the paper is 
to understand the output of Monte Carlo calculations of lattice QCD. All results, however, 
are for matrix elements defined at a fixed (Euclidean) time, so they apply equally well to 
the Minkowski theory. Indeed, with the conventions introduced here, the formulae in this 
paper hold for both kinds of time, unless specifically noted. 

The Euclidean action can be written S = — / d x where C is the Lagrangian, and 
the weight factor in the functional integral is then e~ s . The metric is 5^, Greek indices 
run from 1 to 4, and Dirac matrices satisfy {7^,7^} = 25^. A convenient basis is given in 
Ref. |IJ, in particular 
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(3.1) 

As usual, we take real (Minkowski) time to be t = x°. Then Euclidean time x 4 = ix°, and 
the general rule relating the fourth component to the zeroth component of a four-vector q is 

q 4 = iq°. (3.2) 

Because the spatial components are the same, it is convenient to put all modifications into 
the time component. Therefore, this paper uses the metric g^ u = diag(— 1, 1, 1, 1), where 
Greek indices run from to 3, so go — ~Q° and q 2 = —{q ) 2 + q 2 - And the Dirac matrices 
7° = — i'-f 4 and 7 J differ by a factor of —i from those of the most common Minkowski 
convention. 

The four-volume element is defined to be 

d 4 x := dx 1 dx 2 dx 3 dx 4 = i dx°dx 1 dx 2 dx 3 . (3.3) 



The factor of i in Eq. ( |3.3| ) is the most unusual convention introduced here, but it allows 
many formulae given below look the same in both Euclidean and Minkowski space-time. For 
example the weight factor of the path integral is always 

The foregoing conventions can be used in any field theory. In HQET one introduces a 
velocity v, with v 2 = —1. Although heavy-quark symmetry of lattice gauge theory is only 
guaranteed in the rest frame v = 0, it is convenient to keep v arbitrary. The projectors 

P ± (v) = l(lTifl (3.4) 

project onto "upper" and "lower" components of spinors. For any vector q the components 
orthogonal to v, 

ql = q a + v a v ■ q, (3.5) 

play a special role. In the rest frame they are the spatial directions. It is also convenient to 
introduce 

ri a p = 5 a p + v a v f} (3.6) 

to project out orthogonal components of a tensor, e.g., qj_ = f] a ^ ■ 

In HQET heavy quarks are represented by a heavy-quark field satisfying 

hp>=P + (v)hp\ (3.7) 
The anti-quarks are represented by = P_(v)h^~\ As in the usual HQET one can either 



consider the anti-quarks to be decoupled pqjl9[| or integrated out p0| . But in this paper, 
having shown that the heavy-quark symmetries hold in lattice QCD, the effective Lagrangian 
is developed principally on the basis of symmetry. The heavy-quark Lagrangian is written 



-HQET 



£(0) + £(l) +jC (2) + ... ; (3-8) 
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where the leading term is 



= hp\iv ■ D - (3.9) 

A non-zero rest mass mi is introduced to describe the exponential fall-off of Euclidean Green 
functions, e - ^ 4 with energies E m rri\. The further interactions contain operators of 
dimension 4 + s. By dimensional analysis their coefficients, of dimension — s, contain powers 
of the short-distance scales I/tuq or a. 

The Lagrangian is the unique scalar of dimension four satisfying the Isgur-Wise 
symmetries [BJ. The heavy-quark spin symmetry is manifest, but with m\ 7^ the flavor 
symmetry is not. It is, however, there. In Eq. ( |3.9| ) let the field to be a column vector 
for all the flavors of velocity v, and let mi denote a mass matrix. For example, for two 
flavors 



(3.10) 




Let 9 = (mi c — m\^)v ■ x and consider the generators 

^ =0 • ^n ,^ n > r 3 = -L .1 h ( 3 - n ) 




J' 2 V ze~ lt/ 

satisfying the SU(2) algebra [r a! ,r e ] = e'® e T* . Then the flavor symmetry is 

h (+) e^-ftW, U+) ^ /i(+)e- Ttl ^. (3.12) 

The symbol T>^ = Z>' — imiv 1 *, which was introduced in Ref. [^j, satisfies r d ] = and 
is, thus, trivially covariant under the transformation ( |3.12j ). Therefore, flavor-symmetric 
operators take the form 

QW-Mn = h{+) T V^ ■ ■ ■ V^hl + \ (3.13) 

where T = P + (v)TP + (v). Spin-symmetric operators have T = 1 (or tf, since i/)h^ = ih^). 
The only flavor- and spin-symmetric scalar at dimension four is hj^'iv ■ Vh^\ which is C^°\ 
Thus, the symmetries of HQET with non-zero rest masses are the same as without. 

In the following anti-quarks are not considered further, so from now on the heavy quark 
field is written h v instead of h^' . 

To describe deviations from the symmetry limit, one introduces the higher-dimension 
interactions C^ s \ which are built from operators like Or- These are general enough to include 
the gluon field strength, because F^ u = [D^,D U ] = [D^,T) V ]. One may omit operators that 
would vanish by the equations of motion of C^°\ —iv ■ T>h v = 0. Such operators make 
no net contribution on the HQET mass shell, so they do not appear in on-shell matching 
calculations. At dimension five there can be two T>s, so 

+ < 314 > 

where 
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2 = h v D 2 ± h v , (3.15) 
O b = h v s a pB af3 h v , (3.16) 

with s a /3 = —io a p/2 and B a/S = r^r^F^ . In the rest frame, O2 gives the kinetic energy and 
Ob the chromomagnetic interaction. At dimension six, with three Vs, 

£ p) = ;§- + ;§-, (3-17) 



where 



O c -^[DL^a]^, (3-18) 
E = -h v i<T afi {DZ,iE f> }h v , (3.19) 

with E' 3 = -v a F af3 $ In the rest frame, O d gives the Darwin term and Oe the spin-orbit 
interaction. The complete list of dimension-six interactions includes four-quark operators, 
such as q^qlivV^hv, but their coefficients all vanish at the tree level. 

Distinct inverse masses l/m 2 , l/m^, l/m 2 D , and l/m E are introduced as a notation for 
the coefficients of the modified HQET. One could have equally well written zb/iti2 instead 
of 1/rriB, and so on, but to trace the effects of the higher- dimension operators on physical 
observables the notation of inverse masses is adequate. The numerical factors and powers 
of the inverse masses have been chosen so that all masses become the same in the tree-level 
continuum limit. At non-zero lattice spacing and in the presence of radiative corrections, 
this is no longer guaranteed. 

Concrete expressions for the coefficients lie beyond the scope of this paper. They depend 
on couplings of the lattice action, the velocity v, and the HQET renormalization scheme. 
Ideally one would like to devise a non-perturbative scheme for computing the coefficients, 
but so far they have been studied only in perturbation theory. For the lattice actions in 
common use, expressions are available at the tree level for mj, l/m 2 , and 1/m^ jlj, and at 
the one-loop level for mi and 1 / m 2 |2T| . 

Through dimension six the effective heavy-quark Lagrangian is rotationally invariant. 
Starting with dimension seven, this is no longer the case. For example, consider the term 

£(3) = . . . + a s W4 j- h^Dfhi + \ (3.20) 

i=l 

written in the rest frame, v = 0. In the usual HQET, rotational invariance of continuum 
QCD implies u> 4 = 0. With lattice QCD, however, w 4 does not vanish unless the lattice 
action has been improved accordingly. 

To describe electroweak transitions among hadrons containing a single heavy quark, 
HQET introduces effective operators for the interactions mediating the transitions. Even 
in simple cases, such as the vector and axial vector currents examined below, the num- 
ber of operators in the heavy-quark expansion is large, and the details of the construction 
are different for heavy-to-heavy and heavy-light transitions. The notation for currents is 
postponed, therefore, to Sees. |V11| and |V111| . 



3 The chromoelectric field of Ref. |Q] is related (in the rest frame) to the one here by E [jjj = iE. 
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IV. PROPERTIES OF 



The previous two sections establish that the heavy-quark limit of lattice QCD can be 
described by the effective Lagrangian A -*, with small corrections from 

A = C m + A 2) + • • • . (4.1) 

This means that the eigenstates of lattice QCD are not very different from the eigenstates 
of the quantum field theory defined by Aight + £ > where Aight is the (Symanzik effective) 
Lagrangian of the light quarks and gluons. Apart from the rest mass m\ and lattice artifacts 
of Aight, this is the same lowest-order Lagrangian that is used to describe heavy quarks in 
continuum QCD. 

To use HQET to connect lattice QCD to continuum QCD, one must understand how the 
rest mass and the higher- dimension interactions influence observables. This section shows 
that the eigenstates of the Hamiltonian corresponding to Aight + £hqet are independent of 
the rest mass. In particular, the eigenstates of Aight +&°* do not depend on the heavy flavor 
at all. The remainder of the paper then develops perturbation theory in A around these 
flavor-independent states and studies how the perturbations affect several observables. 

To show that the rest mass m\ decouples from non-perturbative observables, it is con- 
venient to switch to the Hamiltonian formalism of HQET. The canonical conjugate to the 
field hv is |3 

7r„ = iv°h v (4.2) 

so at equal times (x° = z°) 

{h v (x),7r v (z)} = {h v (x),tv°h v (z)} = i5^\x - z)P+(v). (4.3) 
The Hamiltonian H = J d 3 x 7i has the density 

n = ftught + ^ (0) - a, (4.4) 

including a term for the light degrees of freedom. The leading heavy-quark Hamiltonian 
density 

H {0) = TT v d h v - A 0) (4.5) 
= mih v h v + iv°h v A°h v — ih v v ■ Dh v . (4-6) 

From Eq. (|4.3|) one can see that / d 3 x h v h v commutes with all other terms in H, including 
with ifught and Lj — j d 3 x£i. Thus, the eigenstates of H are independent of m\. This 
result is well known in other approaches to the heavy-quark limit |IUL[23|], but the general 
proof within HQET does not seem to be widely appreciated.^ 



In specific examples, a small rest mass, called a residual mass, has been shown to drop out of 



the 1/rriQ corrections [22]. 



This result has a very important consequence. In the HQET description of lattice QCD, 
lattice-spacing dependence appears in three places: the rest mass, the short- distance coeffi- 
cients of £/, and the light degrees of freedom. Because the rest mass drops out of physical 
observables, it is acceptable — perhaps even advisable — to tolerate a discrepancy of the rest 
mass from the physical mass. Genuine lattice artifacts of the heavy quark stem from de- 
viations of the /wp/ier-dimension short-distance coefficients from their continuum limit, and 
the couplings of the lattice action should be tuned to minimize them. To make this more 
point concrete, the effects of Cj can be propagated to observables with tools developed for 
the usual HQET, as shown in the rest of this paper. 



V. PERTURBATION THEORY IN £j 

The previous sections have established that lattice gauge theory with heavy quarks can 
be described by the effective Lagrangian Aiqet, whose eigenstates are close to those of 
the leading-order theory with Lagrangian £(°\ To trace the effects of the higher-dimension 
operators in £j on observables, they can be treated as perturbations. A formalism for 
perturbation theory that exploits heavy-quark symmetry is reviewed in this section. 

When proceeding to second order in Cj, as in Sees. |VT| and [Vn| below, one must be careful 
to be consistent, for example about the normalization of states. Thus, the discussion starts 



(Sec. |V AT ) with a careful setup of time-ordered perturbation theory, to generate heavy-quark 
expansions based on the eigenstates of C^°\ These states are desirable not only because they 
form mass-independent multiplets under heavy-quark symmetry, but also because they are 
affected by the lattice only through the light degrees of freedom. The formalism makes 
no explicit reference to the short- distance coefficients of the modified HQET, so it applies 
equally well to the usual HQET and could be used there as well. The heavy-quark expansion 
becomes a series of terms consisting of short-distance coefficients multiplying matrix elements 
of time-ordered products in the eigenstates of C^°\ Many relations among these matrix 
elements follow from heavy-quark symmetry, and Sec. |V B| reviews the trace formalism, a 
technique for deriving such relations. 



A. Time-ordered perturbation theory 



The perturbative series can be generated by generalizing the interaction picture for vac- 
uum expectation values to transition matrix elements. There are three quantum field theories 
to consider: the underlying theory [here lattice QCD with action ( |2.1| )1; the full HQET with 
Lagrangian (|3.8| ); and the leading HQET with Lagrangian (|3.9|) . The states treated here are 
hadrons with one heavy quark. The (lattice) QCD state with a heavy quark of flavor b (c) 
is denoted \B) (\D)). The analogous full HQET state is denoted \B V ), where the subscript 
labels the chosen velocity. Finally, the infinite-mass states are denoted \b v J;ja), where b is 
the heavy flavor in the HQET with velocity v, J is the hadron's spin, j is the spin of the 
light degrees of freedom, and a encompasses all other quantum numbers of the light degrees 
of freedom. By heavy-quark flavor and spin symmetry, the spatial wave functions of these 
states do not depend on b or J. 



10 



By the Gell-Mann-Low theorem [25[] the lowest-lying (i.e., a = 0) spin- J hadron is 
related to the corresponding infinite-mass state by 



\B) = lim 4 /2 ^T^MiQ) (51) 
1 v) r-oo(i-io+) (6 v J;iO|i7(0,TT)|6„J;iO)' l ' J 



where -Zg is a state renormalization factor and 

rt 

'to 



C/(t, t ) = T exp /* d A x d (5.2) 

Jt t) 



is the familiar interaction-picture propagator. The state renormalization factor has the 
usual interpretation of the overlap between the unperturbed and the fully dressed states: 
Z 1 J 2 = (b v J'J0\B v ). 

To derive the heavy-quark expansion without ambiguities stemming from the normal- 
ization of states, one should set up perturbation theory so that Zb does not appear. For 
example, the energy of a fully dressed state can be written p6l 



(b v J;jO\H\B v ) _ (b v J;jO\H(0)U(0,-T)\b v J;jO) 
(b v J;jO\B v ) (b v J;jO\U{0,-T)\b v J;jO) ' 1 ' ' 

in which the normalization of states clearly cancels. Similarly, matrix elements for flavor- 
changing transitions can be expressed 

(c v ,J' ]J 0\TO 1 ---O n eJ dixCl \b v J; J 0) 



(D v , 


\T0 1 ---O n 


B v ) 


(A/I 


D v i) 1 / 2 (B V \B V )V* 



(5.4) 

where the upper (lower) sign of Eq. (|5.1|) is used for the initial (final) state, and products of 

U(ti,t 2 ) have been coalesced according to its well-known properties |[27|| . The factors Z l J 2 
i li 

and are eliminated in favor of the denominators by taking the modulus of each side of 
Eq. ([TQ). 

The operators Oj in Eq. ( |5.4|) are operators of HQET. In general an operator from the 
underlying theory is described by a sum of operators in HQET, cf. Sees. |V11| and |V111| . On 
the left-hand side the operators have the time dependence 

Oj(t) = e iHt O je - iHt (5.5) 

of the Heisenberg picture, whereas on the right-hand side they have the time dependence 

Oj{t) = e iHot O je - lHot (5.6) 

of the interaction picture. They are related by {H) {t) = U Jf ( y t,0)O (I) (t)U( y t,0). When O 
contains explicit time derivatives, as in some cases in Appendix the time dependence of 
the Us generates additional contact terms in the T-product in the interaction picture. 

This setup of time-ordered perturbation theory is equivalent to Rayleigh-Schrodinger 
perturbation theory P5 ]. The denominators on the right-hand side of Eq. Q5.4p are rarely 
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made explicit in the literature on HQET, but they are necessary. Indeed, in tracing the 
equivalence to Rayleigh-Schrodinger perturbation theory, one sees that the denominators 
generate wave-function renormalization and remove \b v J;jO) and \c v rJ';jO) from sums over 
intermediate states. The procedure is analogous to taking connected vacuum correlation 
functions.^ The infinitesimal in the limit T — ► oo(l — i0 + ) is needed to dampen the integrals 
in Minkowski space, and it is unnecessary in Euclidean space. There is no issue of ana- 
lytic continuation here: the symbol U(t,to) is just an integral representation of the energy 
denominators in ordinary perturbation theory. 

The principal advantage of Eqs. Q5.3| ) and ( p.4|) is that they separate cleanly how each 
term in the heavy-quark expansion affects the matrix element on the left-hand side. As 
desired, the normalization conditions on the full and infinite-mass states cancel separately. 
Each operator in Ci can be treated one insertion after another, and the expansion leads 
to matrix elements in the mass-independent states \b v J;jO) and (c v rJ';jO\. On the other 
hand, a formalism that starts with vacuum expectation values of time-ordered products and 
proceeds to the left-hand side via the reduction formula leads to expressions with "in" and 
"out" states whose masses equal those of the fully dressed states. 

When employing HQET to describe lattice QCD it is especially helpful to obtain a series 
in the mass-independent eigenstates of + £n g ht- These states depend only mildly on 
the lattice spacing, through the light degrees of freedom. Thus, the discretization effects of 
the heavy quark are truly encapsulated into the short- distance coefficients of £/, and one 
can estimate their effect simply by comparing the heavy-quark expansions of continuum and 
lattice QCD. With the expansions derived in subsequent sections, the comparison is made 
easily by substituting the usual coefficients for the modified ones. 

Although not strictly necessary, it is convenient to choose normalization conditions for 
the states. In the underlying theory we normalize plane-wave states so that 

(B(p')\B(p))=v°(2nf5(p'-p), (5.7) 

where v° = yl + v 2 . In continuum QCD v = p/M is the physical velocity of the true 
hadron, and in lattice QCD the relation between v and p should tend to the same as 
pa — > 0. Equation ( |5.7| ) is convenient because it is relativistically invariant and its infinite 
mass limit is well behaved. We also normalize full HQET states so that 

(B v (k')\B v (k)) = v°(2n) 3 5(k' - k), (5.8) 

where k^ is a small residual momentum ||23|| , and likewise for the infinite- mass states. Note 
that in Eq. (|5.8|) the factor of v° does not introduce mass dependence; in HQET the velocity 
is an ingredient in the construction of the effective Lagrangian, not a property of the states. 

To regulate 5 functions one should smear plane-wave states into wave packets before 
expanding out Eq. ( |5.3| ) or ( |5.4| ). With the same normalization condition ( |5.8| ) for fully 
dressed and infinite-mass HQET states, the factors of v and the smearing functions cancel 
completely after expanding. One can thus re-write Eq. Q5.4| ) 



Reference notes both the significance of the subtractions and the analogy with connected 
vacuum amplitudes, but prefers not to use HQET. 
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{D v ,\TO x ---O n \B v ) = (c v ,J'-j0\TO 1 ---O n ef d4xCl \b v J; 3 0y (5.9) 

where the star on the right-hand side is a reminder to include the extra terms generated by 
expanding out the denominator of Eq. (|5.4T) . In Sec. |VI]] and Appendix |A| this notation is used 
for T-products (<v J'\T OO b x \b v J}* , (c v > J'\T O x OO b Y \b v J)\ etc., where the operators O h x are 
those appearing in Lj for flavor h. The star means to collect all terms from the expansion 
with the specified insertions. 



B. Trace formalism 



To evaluate the right-hand side of Eq. ( |5.9|) there is a powerful formalism, called the 
trace formalism, which takes full advantage of heavy-quark symmetry P9"| . The objective is 
to calculate transition amplitudes of the form 

T h A Xc A f = (c v 'J , ;MTh v/ T 1 Gt 1 h v/ ---h v ,T n G^h v ---h v T N Gi f N h v \b v J ] jOy (5.10) 

and 

T b % A * = (0\TqT l G^h v h v T 2 G^h v ---h v T N G^h v \b v J;jO)\ (5.11) 

where the G Ak is a combination of covariant derivatives T> (including field strengths F^ u ) 
and light-quark bilinears qq with Lorentz indices abbreviated by the superscript A^. 

The color and spin dependence of each static propagator Th v (x)h v (y) [or Th v >(x)h v '(y)) 
factors into a Wilson line and a projector P+{v) =: P + [or P+(v') =: P' + }. That means that 
the amplitudes can be written (for j = \ mesons) 

T h %y = - ti{MAv r )T 1 P[ - ■ ■ P[T n P + - ■ ■ P + T N Mj(v)E A ^}, (5.12) 

and 

r b % AN = -ti{Y l P + Y 2 P + ■ ■■P + r N Mj(v)~ A ^ AN }, (5.13) 

where Aij and fAj* are spin wave functions and H Ar " An parametrizes the spatial wave 
functions and a trace over color of the Wilson lines, punctuated by the G Ak , with the light 
quark propagator. There is only one trace over heavy-quark spin, because products of traces 
correspond to disconnected terms, which are subtracted when expanding Eq. ( |5.9| ). The 
minus sign arises because the trace over spin is obtained after anti-commuting the left-most 
quark field all the way to the right. 

Spin wave functions such as M.j{v) and A4j>(v') are determined by spin symmetry alone. 
For j = \ they are 

M (v) = t2- 1 / 2 P + (v) l5 , (5.14) 
Mi(v) = i2~ 1/2 P + (y)4. (5.15) 

Charge conjugates are M. = 74^^74, 
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M (v) = t2- 1 / 2 l5 P + (v), (5.16) 
M 1 {v) = i2- 1 ' 2 4P + {v) ) (5.17) 

where e = e* in Minkowski space-time and e = (e*, —64) in Euclidean space-time. Note that 
M = P+M.P- and M = P-MP+. Generalizations to j — and j = 1 baryons [050 anc ^ 
to higher angular momentum El are available in the literature. 



The functions E Al ' " An cannot be obtained from symmetry considerations alone. They 
depend on the velocities v' and v and the quantum numbers of the light degrees of freedom. 
They parametrize the long-distance dynamics of C^ + Cn g ht, so they do not depend on flavor, 
and they suffer from lattice artifacts only of the light degrees of freedom. As explained above, 
cutoff effects of the heavy quark are captured in the coefficients of the modified HQET, which 
multiply matrix elements (|5.10| ) and (|5.11|) . 



VI. HADRON MASSES 

The simplest application of the HQET formalism is to generate an expansion for the 
rest mass of a heavy-light hadron. In numerical lattice calculations the energy of a state 
of momentum p is computed by looking at the (imaginary) time evolution of a correlation 
function 



(M*4)*i(o)> = Vp 



^(x 4 )^e-^"^Ki?J^|0)| 2 + ^-X4)E^' (p) l(5n'|$ P |0)p 

n n' 

(6.1) 



where \B n ) i\B n ')) are full lattice-QCD states connected to the vacuum by (^p)- By 
a combination of judicious choices of $^ and taking 24 large enough, one can isolate the 
lower-lying states. At small momentum, the relation between energy and momentum is 

V 2 

E(p) = Mi + — , (6.2) 
Kyj 2M 2 ' y J 

which defines the hadron's rest mass Mi and kinetic mass M 2 . (Some authors call Mi the 
"pole" mass, but Mi and M 2 are both properties of the particle's pole.) In this paper 
upper-case is used to denote hadron masses, and lower-case to denote quark masses. 

These energies can be thought of as eigenvalues of a Hamiltonian, defined via the transfer 
matrix, which HQET models with Eq. (|4.4j ). In Eq. ( |5.3|) H is always to the left of U(0, —T), 
so one can make the split H = Hn^t + — Lj and act the first two terms on the 
bra (b v J;jO\. Setting p = and calling the leading eigenvalue 

, a _ (KJ;jO\[H nght + H^}\b v J;jO) 

the heavy-quark expansion of the hadron mass is generated by 

Mi = mi + A - (bvJ-JOlLjTef^lbvJ-jay, (6.4) 
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where Lj is at time 0, the time integration is from — oo to 0, and the star is a reminder 
not to neglect the denominator in Eq. ( |5.3|) . The quark's rest mass enters solely additively 
because its term in the Hamiltonian commutes with all others. 

The expansion of Eq. ( |6.4|) leads to reduced matrix elements that depend on the spin j 
of the light degrees of freedom (j = for the Af, baryons, j = 1/2 for the B and B* mesons, 
etc.), but not on the heavy quark's spin. Through order 1/mg one defines 



(b v J;j0\O 2 \b v J;j0) = \ u (6.5) 

(b v J;j0\O B \b v J;j0} = dj\ 2 , (6.6) 

(b v J;j0\O D \b v J;j0) = -2p u (6.7) 

(b v J;j0\O E \b v J;j0) = -2d jP2) (6.8) 

and, in the notation of Ref . J^J , 

J d 4 x(b v J;j0\O 2 (0)O 2 (x)\b v J;j0y = T x , (6.9) 

J d 4 x (b v J;j0\O B (0)O B (x)\b v J;j0y = % + dj(% - T 2 ), (6.10) 



d 4 x (b v J;j0\O 2 (0)O B (x)\b v J;j0y = / d 4 x (b v J;j0\O B (0)O 2 (x)\b v J;j0)* = djT 2 . (6.11) 



The J-dependence in Eqs. ( |6.5|) - (|6.11|) is do = 3 (for the B meson) and d\ — — 1 (for the 
B* meson). For the baryon there are fewer non- vanishing matrix elements; the above 
formulae hold if one sets di/ 2 = 0. The parameters A, A n , p n , and T n are the same as in 
continuum QCD, apart from lattice artifacts of the light degrees of freedom. Combining 
Eqs. (|6.4j )- (|6.11| ) the rest mass becomes 

Ai dj\ 2 pi djp 2 Ti 2djT 2 % + dj(% - T 2 ) 



2m 2 2m B Am D 4m E 4m 2 2m 2 2m B Am B 

(6.12) 

The result ( |6.12|) is simple enough that it could have been written down upon inspection of 
Eqs. ( |3.14j ) and ( |3.17| ) and comparing to the continuum papers p3|j3^j28|j26|1 . 



This result is the first example of the expansion for which Eq. (|1 . 1| ) is a prototype. 
Short- distance effects of the heavy quark, including lattice-spacing effects, are contained in 
the "masses" m\, m 2 , m B , m^, and m B - If the bare mass is adjusted so that m 2 = mq, then 
the mass formula ( |6.12| ) shows that the spin- averaged splittings, such as m^ b — ^(m B +3m B *), 
are reproduced correctly to order 1/rriQ. The Sheikholeslami-Wohlert action has a second 
parameter, with which l/m B can be adjusted (via a short-distance calculation) to reproduce 
correctly the spin splittings, such as m B * — m B , to order I /nig. These adjustments are 
essential, because in matrix elements the rest mass plays no role whatsoever. 

In the usual HQET with mi = 0, the quark mass is added to A and the higher-order 
terms. Ambiguities of the HQET renormalization scheme, including those of infrared renor- 
malons in the on-shell scheme, cancel in the sum. Similarly, the difference m 2 — mi can be 
added to Eq. ( |6.12|) : M = Mi + m 2 — mi. Adding the residual mass in this way has the 



virtue that m 2 — mi does not suffer from infrared ambiguities, even in the on-shell scheme. 
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VII. SEMILEPTONIC FORM FACTORS 



Another interesting application of HQET is the heavy-quark expansion of form factors 
in the exclusive semileptonic decays B — > D*lv and B — > Dlv. These decays offer the most 
promising way to decrease the uncertainty in the CKM element \V c b\, provided the hadronic 
matrix elements can be calculated reliably. Recent work [^,^] shows that calculations of 
the form factors at zero recoil with statistical errors of a few percent are feasible. The aim 
of this section is to describe the l/mq and I/ttlq contributions to the lattice observables 
calculated in Refs. |]35|j36| , and compare them to the description of the form factors in the 



usual HQET. The technical details are in Appendix [A], mostly following Refs. |3"5| , |3~3 
The transitions are mediated by the charged weak currents 

V = ci^b, A" = c$7"75&, (7.1) 

where c and b are conventionally normalized continuum quark fields. Currents in lattice 
gauge theory and in HQET are introduced below, but the symbols V^ 1 and A! 1 are reserved 
for the physical currents. The hadronic part of the transitions involves the matrix elements 
(D(*)|V*|S> and (D*\A^\B). For B -> Dlv there are two form factors h + and h_. With the 
normalization ( |5.8|) they are related to the matrix element by 



(D(v')\V»\B(v)) = \{v' + v)»h+(w) - \{v' - vYh-iw), (7.2) 

where w = —v' ■ v. Zero recoil corresponds to w — 1. In Eq. ( |7.2| ) the final velocity is kept 
distinct from the initial velocity to be able to obtain h-(l). For B — > D*lv there are three 
axial form factors, defined by 

(D*(v'^)\A»\B(v)) = \{w + l)iTh Al {w) + li7-vv^h A2 (w) + li7-vv ,tM h A3 (w), (7.3) 

and a vector form factor, but at zero recoil the decay rate depends only on h Al (l). For 
reasons that will become clear below, the zero-recoil matrix element 

(D*(v,e')\V»\B*(v,e)) =7-ev> l h 1 (l) (7.4) 

and its form factor hx(l) are also of interest. 

Note that continuum QCD currents define the form factors. To generate the heavy-quark 
expansion of these form factors, one replaces the currents V M and A^ with effective currents 
built from the heavy-quark fields and the fields of the light degrees of freedom. The effective 
currents and the heavy-quark Lagrangian are treated to the desired order in 1/rriQ, and 
Eq. (|5.4[) should be used to generate the expansion, consistent to that order. 

The zeroth order is simple and worth reviewing briefly. The QCD currents are related 
to HQET currents via 

V M = rjvCyd^bv - IflvW - vYc v >b v - ^ v (v' - v) u c v da^b V) (7.5) 
A 1 * = VACv'h^^bv - \Pa{v' - vYcy'^bv ~ \ia{v' - v) u c v >ia^-f 5 b v , (7.6) 

where the symbol = means that the operators, though defined in different field theories, have 
the same matrix elements. The short-distance coefficients depend on the two masses; rjj and 
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7j are symmetric upon interchanging the masses (j G {V, A}); (3j is anti-symmetric; at the 
tree level they satisfy t]j = 1, (3j = jj = 0. To obtain the leading heavy-quark expansion, 
one simply takes matrix elements of the effective currents in the states of the infinite-mass 
theory. From the trace formalism one finds 



h+(w) = [rj v + \{w - 1)tv] £(tu) + 0(l/m Q ), 
h-(w) = \{w + l)/3v£(w) + 0(l/m Q ), 
hi(w) = [r] V + \{w - 1)tv] £(to) + 0(l/m Q ), 
h Al (w) = T)a€(w) + 0(l/m Q ), 



(7.7) 
(7.8) 
(7.9) 
(7.10) 



with a single HQET form factor called the Isgur-Wise function. At zero recoil it 

is normalized by heavy-quark symmetry ||, so £(1) = 1. Therefore, the leading term in 
heavy-quark expansion is h + (l) = h\(l) = r] V , h_(l) = f3 v , and ^^(l) = f]A- 

The 1/rriQ p9| and 1/m.Q |3^j3^| corrections to Eqs. (|7.7|) - (|7.10|) have been worked out 
with HQET. This section repeats the analysis through order I/tuq for the lattice approxi- 
mants to the form factors introduced in Refs. p5| , p6| . The only crucial difference is that the 
short-distance coefficients are tracked carefully and their contributions are kept separate in 
the final results. 



A. Lattice and HQET currents 



To compute the form factors in Eqs. (|7.2| )-(7Ti) with lattice gauge theory one introduces 



combinations of lattice fields with the same quantum numbers as V M and A^. The lattice 
currents are given by a series of dimension-three, -four, -five, etc., operators, with coefficients 
chosen to attain the right normalization and to reduce lattice artifacts. Several choices have 
been made in the literature, but with Wilson fermions they can all be described by HQET: 
the different choices simply have different short-distance coefficients. 

Let ZycbV^ (Z^cbA^) denote the lattice approximant to the charged b —> c vector (axial- 
vector) current. To conform with much of the literature on lattice gauge theory, the current's 
normalization factor Zjcb in shown explicitly. Then, suppressing the space-time index, the 
lattice currents are related to HQET currents via 

z V c b v lat = + E E v(r,s ~ r) ( 7 - n ) 

8=1 r=0 

= V (P) + 00,1) + ^(1,0) + y(0,2) + y(l,l) + y(2,0) + . . . ^ ^ 

and similarly for A(^ t . The HQET operator V^'^ carries dimension 3 + r + s. To make 
contact with the usual HQET, it is helpful to think of the dimensions being balanced by 
r powers of l/m c and s powers of The dimension-three vector current is 

V; (0) = (W + Sv l ^)c v n lfl b v - - v)^ v ,b v - ^Jv' - vYc v na, u b v . (7.13) 

In general the coefficients depend on the directional indices, because the lattice singles out 
the time direction. The overall factor Zy is conventionally chosen so that 5rjy^ = 0. Then 
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5r]y^ vanishes at the tree level, but not in general. As with the usual HQET /3y and 7^ 
are, respectively, antisymmetric and symmetric upon interchange of heavy quark masses and 
both vanish at the tree level. The operator multiplying /3y (7^) makes a contribution at 
first (second) order in v ' — v. 

At dimension four and higher many operators arise, and a complete catalog requires a 
voluminous notation. Only the 77-like terms are listed here. /3-like terms are not needed until 
Sec. |V11D| , and 7-like terms are not needed at all. With this restriction, the dimension-four 
currents are 

(714) 



T/(1,0) _ 1 „(1,0) Cy'P^^K 
v it — ~r f V ^ • 



Vv 2m 3c ' (7 - 15) 

where D±> = D + v ' ■ Dv'. The notation t]y' ^ /m 3c and rjy' 1 ^ /m^ for the short- distance 
coefficients follows a helpful convention: for degenerate quarks the coefficient is merely 
l/m.3, which thus depends only on the indicated flavor; rjy then describes the additional 
radiative corrections for non-degenerate masses. The dimension-five currents are 

t/(0 2)_ (0,2) Cyn^Dlby t (o,2) c v >i^^s al3 B a pb v (0,2) c v >i'y fJ ,ifib v (>? -ia\ 

V -^VDl o m 2 +VVsB 2 VVaE 7 2 ' ^' A ' 3 ) 

° m D\b ° m sBb ^ m aEb 

T/(2,0) _ _(2,0) ^Dj^TjA; (2,0) ^yS^B^^A (2,0) Cy^^A (rj , ? , 

± Xm D 2 c om sBc 4m aBc 



y(i,i) _ (P±'il^P±)iby _ c v ' (P±'i / y fl p A _)sb v 

M yi 2m 3c 2m 3b Vs 2m 3c 2m 3b 

(r,s) 



(7.18) 



where again 1 / m 2 ^ depends only on the indicated flavor and r/y depends on both masses. 
The two coefficients Zy^ and z^ 1 -* multiply the spin-independent and spin-dependent part 
of the Dirac matrix structure. They do not reduce to 1 for equal masses, because l/m 3 is 
defined through the dimension-four currents, but for most choices of the lattice current they 
do equal 1 at the tree level. 



B. At zero recoil: /i+(l) and h\(l) 

The matrix elements that are to be described are 

(D\Z jJlat \B) = (A/b' (0) |A,> + (D v ,\jM\B v ) + (D v ,\jW\B v ), (7.19) 

where j is V or A, and = j (0il) + j (1 ' 0) , j {2) = j (0,2) + + j (2,0) . The first two matrix 
elements on the right-hand side of Eq. ([7.19Q must be expanded via Eq. (|5.4|) to second and 
first order in £/, respectively. There are, consequently, many HQET matrix elements to 
introduce. The matrix elements and their abbreviations, analogous to those in Sec. [VT] , are 
listed in Table |. The notation mostly follows previous work [0,0 . 
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TABLE I. Notation for HQET matrix elements in Refs. [33, 34 and this work. 



contribution 


Ref. pi 


Ref. pi 

1 1 


this work 


<j (0) > 




1 


£(w) 


0' (1) ) 


^±,3(^) 








Ai(w) 


Xi 


Ai(w) 


(Tj®0 B )* 


^2,3^) 


X3 


A 2j s(w) 


(J' (2) ) 


</ ) 0,...,3(u ; ) 


Al, 2 


Ai,...,4<y) 








S 2 3 (u>); -Fi i(w) 


(TjWO B )* 


£?4 ii (w); £vi ii (u>) 

j A J. v / " 41,..., IX \ / 




S4,...,ll(u>); i ? 4,...,ll(w) 




2Bx(w) 


2Ei 


Ai 




2^2,3^) 


2H 3 


A 2 


(T0 2 jMo 2 y 


£>iH 




L> 


(T0 2 jWo B y 


-D 2 ,3(w) 




E 


{TObj^Ob)* 


At,...,ioH 


#1,2 


Rl,2 


(TjWo 2 o 2 y 






A = -\D 


(TjW0 2 B )* 


C 2 , 3 (w) 




B = -E 


(TjWO B B )* 


C4,...,12(^) 


Cl,3 


Cl,3 = — 2"-Rl,2 



One can work out the matrix elements using the trace formalism. At zero recoil 
(D v \j^\B v ) vanishes. For the vector-current transitions B — ► D and B* — > 7J* with 



v = v 



one finds 



(D«|z yc ^|5 w > = ww4 J + ^, 



,(2) 
JJ J 



in which (D V \V<®\B V ) yields 



(0) 

.1.1 



1 + 

+ 
+ 



4m2 C 4m 2 



/ 1 



1 1 

+ 



1 



2b, 



V 2m 2c 2m Bc 2m 2b 2m Bb 
1 1 



djB 



Am Bh ) 2m 2c 2m 2b 

1 1 

+ 



2ms c 2m 2b 2m 2c 2m^ 



djE + 



D 
1 



2m Bc 2m Bb 



and (A,|"^ (2) |A,) yields 



W 



(2) 

,y,/ 



/ (2,0) (0,2) 
^VD 2 ^VD 2 



,(1,1) 
-VI 



V 



5WI 



£> 2 . 



2m 3c 2m 



Ai + 



„(2,0) (0,2) 

'l y sB i 'lvs B 

9 ~T~ ~ o 



3/; 



5m 



sBc 



5WI 



sBfe 



+ djR 2 ] 



2m 3c 2m 3b 



(7.20) 



(7.21) 



rfjA 2 . 



(7.22) 



The subscript JJ' denotes the initial and final spins, although here J' = J. The spin factor 
do = 3 for B ^ D and cfi = — 1 for Z?* — > D*. The coefficient factors reveal the origin of 
the contribution. By heavy-quark symmetry Ai and A 2 are exactly the same as in Sec. [VI], 
and A, B, Ci, C 3 , D, E, Ri, and R 2 are new constants parametrizing the light degrees of 
freedom, introduced in Appendix |A| and the last six rows of Table |. 
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Equation ( |7.20|) gives lattice approximants to the form factors h + (l) and hi(l). One 
striking feature of Eqs. ( 7.20| )- (|7.22j) is that there are no contributions of order 1/itlq. For 

Matrix elements of j^ 1 ' in the 
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continuum QCD, this is known as Luke's theorem 
infinite-mass states contribute only when v' ^ v , so a single power of l/m 3 does not appear 
in Eq. ( [7.20| ) . As shown in Appendix [A 1 b| , terms with a single power of 1/777,2 and I /tub 
are absent as a consequence of heavy-quark symmetry and Eq. ( |5.4|) . Thus, Luke's theorem 
holds for lattice QCD also. 

At order l/rriQ, matrix elements that might have multiplied l/m 2 D and l/m% also vanish; 
so do matrix elements involving four-quark operators. Furthermore, the parameters A, B, 
Ci, and C3 can be eliminated, as indicated in the right-most column of the last three rows of 
Table |. As shown in Appendix |A 1 ej this is another consequence of heavy-quark symmetry 
and Eq. ( |5.4| ). Taking these relations into account 



W 



(o) 
,/,/ 



I A 2 D 



A 2 A B d J E-lA 2 B [R 1 + djR 2 



(7.23) 



where 



2m 



Xc 



2m X b 



(7.24) 



So Wjj is correctly reproduced if I/7772 and 1/itlb are adjusted to their continuum values, 
in particular if the analysis identifies m 2 with the heavy quark mass. 

Another lattice approximant to h + (l) and hi(l) is given by double ratios introduced in 
Refs. E 



Ri 



(D\K.\B)(B\V^\D) 
(DKJ^iBKJB) 



(D*Kt 


B*)(B*\V° t 


D*) 




D*)(B* 


V° 

Hat 


B*) 



are approximated by py Q R +) x, where p Vo 



Then \h +> i{. 

advantage of the double ratios, let us rewrite W 



(2) 
.1,1 



ZycbZybc /Zyc 

(2) 
J J > 



w 



(2) 



,/,/ 



(1,1) 



5W 



(2) 



./,/ 



iA 2 

2^3 
/ „(2,0) 



^1 + Zy'Pdj\ 2 



(7.25) 
(7.26) 
To see the 

(7.27) 



+ 



+ 



(2,0) 
VVsB 



+ 



-.(0,2) 

8m 2 Dli 

(0,2) 
VVsB 



(1,1) 
Z V1 

8m| c 
,(i,i) 
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From Eq. ( |7.20|) and the definitions one finds 

PVo jR^ l = r ]v wf}+wf} 



(7.29) 



(2) 

The contribution of 5W)j , which stems from the dimension-five currents, largely cancels. 
Hence, the double ratios depend most strongly on l/m 2) 1/tub, and 1/7773, namely the 
coefficients in C^ 1 ' and V,^ 1 '. 



20 



Equation ( |7.29 ) is an important practical result. If one tolerates errors of order a s jm 2 Q 

and 1/wiq, then pv ^R+,i only requires m 2 = tub = m 3 and Zy' 1 ' = 1 at the tree level, 

and details of the currents V^ ' 2 ^ and V^ 2 ' ^ do not matter at all. With the widely used 
Sheikholeslami-Wohlert action 0], this accuracy is easy to arrange Pl,|3"5|. In practice an 



error comes also from r\y and pv , which are available only to two loops |37j and one loop p8 



respectively. So the recent result f35j for h+(l) has a heavy-quark discretization effect of 
order a 2 , which could be reduced by calculating py Q to two loops. 



C. At zero recoil: ^Ai(l) 

To obtain lattice approximants to h Al (l) one must work out Eq. ( |7.19| ) for a B — > D* 
transition mediated by the axial current. In HQET the currents are as in Eqs. ([7.13 )- 
( |7.1S| ) with a factor 75 inserted in the obvious places. In this case, the overall factor Z A cb is 
conventionally chosen so that 5rf^ cb = 0. 

i 

A useful matrix element has the D* spin is aligned along the i direction and v — v' — 0. 
One finds 



(D*\Z Acb A\jB) = r) A *W$ + W$ + 5W{ 2) 



01 



(7.30) 



in which (D*\A^\B V ) yields — after eliminating A, B, Ci, and C 3 



W, 



(0) 



01 



1 - ±A 2 (A 2 D - 2Q B E) - ±A B (A B i2i - Q B R 2 ) 



2m Bc 2m Bb 



(±R 1 + 2R 2 ) (7.31) 



with 







x 



2m 



+ 



Xc 



2m 



Xb 



(7.32) 



As before, the zero-recoil matrix element does not depend on the dimension-six Lagrangian. 
The matrix element (D*\A^\B V ) of the dimension-five current yields 



W, 



(2) 



01 



5W f 



(2) 



01 



IA 2 + i - 

2 3 3 2m 3c 2m 3b 

( (2,0) (0,2) 

— H - 

„(2,0) 
'lA cb sB 



~(l>l) \ 

A cb l 1 



|A 3 e 3 



1 



,(1,1) 

'A cb l 



y (l,l) 



6m 2 2 



D\b 



6m 



3c 



6m 
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o„(0,2) 

OII A cb sB 



6m 



sBc 



6m 



sBb 



2m 3c 2m 3b 



Ax 




(1,1) \ 

Z A^ S A 2, 



(7.33) 
(7.34) 



after grouping terms as in Eqs. ( [7.27] ) and ( |7.28| ). 
Reference |36| introduces a third double ratio 

Ra 1 



(D*\Al t 


B)(B*\Al t 


D) 


(D*\Al t 


D)(B* 


^lat 


B) 



(7.35) 
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After substituting for each matrix element the foregoing expressions one finds 

PaJrZ = f, Acb w£ ] + Wg ] + 0((v { r } - l)/ml), (7.36) 



where p\ = Z^cbZ^bc/Z^ccZ^bb, f] 2 Acb = r\ A cb7] A bc/rj A^ c "i]A bb i and 



W$ ] = 1 - \A\D - A 2 A B E + lAKR, + 3R 2 ) (7.37) 

where ZA cb — T] A cb Z A cb 

/r] A cb. As before the contribution SWqx of the dimension-five currents 
largely cancels and the double ratio depends most strongly on l/m 2 , l/m B , and l/m 3 , 
namely the coefficients in and A^> . 

Note, however, that pa\J Ra x does not yield hAxiX) but 

^r D *(i)^r s *(i) 



Nevertheless, if the action and currents are tuned so that l/m 2 , l/rug, l/m 3 , and 
match the usual HQET (to a desired accuracy), the three double ratios R + , R\, and Ra 1 
can be combined to yield the l/m^ contribution to /^(l). For example, if one tolerates 
errors of order a s /mQ, as well as I/itlq, one only requires m 2 = mg = at the tree level, 

and one may set Zy' ^ — z A ' ^ — 1- Then, dropping the distinction between Tfi2 1 Trig, and 
m 3 , the double ratios are 

p v R + = 1 - 2A 2 £ P , (7.39) 
p v R 1 = l - 2A 2 £ V , (7.40) 
p 2 A R Al = 1 - A 2 (£ P + £ v + £ A ), (7.41) 

where p v = p v /w, Pa = Pa/va, and 

2£ P = D + R 1 - Ai + 3(2E + R 2 - A 2 ) (7.42) 
2£ v = D + R x - Ai - (2E + R 2 - A 2 ) (7.43) 
^ = |(A 1 - J R 1 ) + 2(A 2 -i? 2 ). (7.44) 

In the approximation being considered the desired form factor is 

Wl) =1 _ A (i^_ (7.45) 
77,4 \2m c 2rrib / 2m c 2m\, 

By fitting Eqs. (|7.39|) - ([7.41|) one can extract £p, and £p + £y + and then one has the 
information necessary to reconstitute /^(l). As with h + (l) there are, in practice, further 
errors because tja and pa are available only at finite-loop order [0,^ . 
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D. Near zero recoil: 



To extract h-(l) matrix elements with non-zero velocity transfer are needed, and some 
new features appear in the analysis. For example, lattice approximants to h- receive a 
contribution from the term in proportional to (3yK Similar "/3-like" terms omitted 
from Eqs. ( [7.14|) -( [7.18| ) also make contributions. We shall not write out all these terms but 
indicate instead how they contribute to matrix elements. 

Suppose one extracts the form factor from a matrix element with v' = —v, pointing in 
the i direction. Then 



(D{-v)\Vi\B{v)) 



lat w (0) , v (2) y(l) 
Vi vv 00 t -'00 "r ^-00 



(7.46) 



where Wqq is given in Eq. ( |7.23|) , and Yqq is like Wqq but with /3-like coefficients replacing 



(2,0) (0,2) 



Vv ' i Vv ' i an d z v'^ ■ The expression in braces is a lattice approximant to h-(w). For 
infinitesimal i>, the matrix element {D V >\V^\B V ) yields 



'(2) 



(2) 



X, 



(1) 



4^ 

2m 



(o,i) ' 
Vv 



3c 



2m 
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(1,0) (0,1) ' 

Vv + Vv 



2m 



3c 



2m 3b 



2&(1) - A + 2E 2 H 3 (l) + 2E B S_(i; 
A 2 0o(l) + A B 0_(1)1 , 



(7.47) 



with A 2 , A B as in Eq. (f7T2l]) and 



1 1 

+ 



2m Xc 2m X b 
The chromoelectric part of (D V '\V^ 2 '\B V ) yields 



(7.48) 



A, 



(2) 
00 



„(2,0) 
'iVaE 



4m 2 aEc 4m aEb 



77 (0 ' 2) \ 2 



(7.49) 



The coefficient factors, together with Table |, make clear the origin of each term. The 
infinite-mass matrix elements A, Ai, and A 2 are exactly those introduced earlier, and the 
new ones £3(1), S 3 (l), H_(l), (fio(l), and 0-(l) are introduced in Appendix |A 2| . As before, 
the dimension-six effective Lagrangian drops out, but the dimension-five currents contribute 



in several places: in Wqq\ Yqq\ and Xqq . Further operators 
and Cv'E'^bv, whose coefficients vanish at the tree level, modify the short-distance coefficients 
of A in Xqq 1 and of Ai in X^ . Thus, many short-distance coefficients influence the accuracy 



(2) 



-ic v 'Dj_b V} ic v 'Dj_,b V} c v rE^b v , 



of Eq. flfgp . 

The main drawback of Eq. ( |7.46| ) is, however, the requirement v' = — v for hadrons of 
unequal mass. Numerical calculations employ a finite volume and, hence, discrete momen- 
tum. Moreover, with the many "masses" the relation between momentum and velocity is 
not plain. To remove these ambiguities Ref. [BSJ introduced another double ratio 



(D\V° t \B) (D\V^\DY 



(7.50) 
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In the spatial matrix elements, the initial state is at rest and the final state has a small 
velocity in the i direction; in the temporal matrix elements, initial and final states both are 
at rest. In continuum QCD, the analogous first ratio is [cf. Eq. (|7.2j )] 



(D\V\B) 



2 u i 



Mi) 
h + (i) 



(7.51) 



2 u n 



(7.52) 



(D\W\B) 

to first order in v' i: and the second is 

(D\V°\D) 

because in the elastic case h + (l) = 1 and h-(w) = 0. Thus, with a suitable adjustment of 
the lattice currents, one can use R_ to obtain a lattice approximant to h_(l)/h + (l). 

In the double ratio of lattice currents the (mass-dependent) factors Zy cancel. With the 
results of Appendix |A2| , and noting that rjycc = 1 and f3y% = 0, 



R. 



rjycb + 5r]yi b 



rolat , V V) , y-(i) , v-(2)i 



77^(1 + Srftfec + 5rfy 

i 



+ 0(l/m; 



q)' 



(7.53) 



Here iy ( 2) ,X^, and r ( 2) 



where 5^y c \ is a combination of Srjy S) /r]y' s> and (tfr/y* — Pyf)/rj 

(2) 

are precisely as above, though in the denominator Wq is evaluated with flavor c in both 
final and initial states. As with the other double ratios, one would like to extract the long- 
distance information from To do so one must have a way to calculate the short- distance 
coefficients, either to adjust them so (3y = (3y and Sr^ = 5fjy^ = 0, or to constrain a fit. 
A simple version of the latter strategy is available if one tolerates errors in h_(l) of order 
1 1 ;+ :•• 4 — — rriB at the tree level, one may 



(2) - - 

and one may neglect Yqo arL( l "V 



a s /rriQ and a s /rriQ. Then it is enough to adjust m 3 

Set 7]y' ^ = 7]y'^ = 1 

hand, Eq. (|7.53| ) can be rearranged to yield 

r] V c b [l - (1 + Sift^R-] + 8rfy% - (3- 



m 2 



«( 2 )w( 2 ) 



V 



Wq . In the approximation at 



lat 

ycb 



-^00 + %t 



(2) 
00 • 



Setting m 3 = m 2 = m B , but keeping m a E distinct, Eqs. ( |7.47| ) and. ( |7.49| ) yield 



X$ + ^oo j = + + ~Aa E X aE [\x + 3 A 



-(2) 



(2) 



(7.54) 



(7.55) 



? (i) 

8(2) 



where 

26(1) -A, (7.56) 

2S S (1) + 2S_(1) - O (1) - 0-(l). (7.57) 

One may fit the left-hand side of Eq. ([7.54 ) to the right-hand side of Eq. ( |7.55| ) with 1 / m 2 aE 
at the easily obtained tree level. After the fit one may reconstitute h_(l) from 

Mi) 



(3 V + A£ x) + AS £^> + | (Aa + 3A 2 ) 



h+(l) 



(2) 



(7.58) 



In practice, there are also errors of order a™ because the coefficients rjy, (3y, Sr/y 1 , and fly 1 
are available only to a finite loop order. Note that the matrix element Ai + 3A 2 appears also 
as the I/rag correction to the pseudoscalar meson mass, cf. Eq. fl6.12|) , so a simultaneous 
fit may turn out to be useful. 
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VIII. LEPTONIC DECAYS 



A straightforward application of the trace formalism gives the first-order heavy-quark 
expansion of the matrix element in leptonic decays. The result for lattice QCD is in Ref. |[39|| , 
but for completeness the derivation is given here. 

With the states normalized as in Eq. (|5.8|), the QCD amplitudes appearing in leptonic 
decays of heavy-light pseudoscalar and vector mesons can be written 

(0\A>*\H{v)} = w^ H /V2, (8.1) 
(0\V»\H*(v,e)) = e^ H */V2, (8.2) 

where V M and A^ are now the vector and axial vector currents with a light and a heavy 
quark, and H (H*) is the pseudoscalar (vector) meson with heavy flavor h. The relation 
between the parameter 4>h and the conventional pseudoscalar meson decay constant is 



(j) H = fjjyj Mij. (8.3) 

There are several conventions for defining the vector meson decay constant, but only is 
considered here. 

In lattice gauge theory the decay constants are approximated with matrix elements of 
lattice currents Z Vq hV qh and Z Aq hA qh with the same quantum numbers as V and A 11 . As 
before, they are not made explicit, to allow for a variety of choices. The underlying currents 
are described by HQET currents, 

Z vq hV* h = rjyqhqijuhv + Cv^v^qK - -^-qi>y^p ± h v H (8.4) 

(o,i) 

71 

Z M hAf = rj Aqh qi^^h v + C A qHV v qia ^ 5 h v - -^-qi-f^ 5 p ± h v H (8.5) 

where q is a light anti-quark field. The coefficient l/2m^ is defined through the degenerate- 
mass heavy-heavy vector current, and r/f' 1 ^ captures the remaining radiative corrections. At 



the tree level ?7j°' 1 ' ) = 1. The coefficients Q vanish at the tree level, and the operators that 
they multiply do not affect </>#(*). Additional dimension-four operators, whose coefficients 
vanish at the tree level, are not written out. 

The static limit is given by the matrix element of the first term of the HQET currents: 

(0|gziyi„|/i„J) = -|0ootr [iT^Mj] 

= iu^oo/V^, (8.6) 

where = 7^75 or 7 M and = t> M or — ze M , for J = or 1. The constant 0oo/2 is introduced 
to parametrize the light degrees of freedom; in the static limit, <$h = <Ph* = 4>oo- As with 



the quantities introduced in Sees. |VI] and [VII] , 0^ differs from its continuum limit, but the 



difference stems only from the light degrees of freedom. 

At order 1/rriQ there are three contributions to </>#(*), from the kinetic and chromomag- 
netic energy, and from the correction to the current. They take the form 
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\Zjjn\HW) = TijiOlqiTMhvJ) + ^J d4x (0 1 T qiV^hv (0) 2 (x) | h v J) 



+ 



Vj 



2m 2 

d A x (0|T qiV^K^OBix^KJY 



•7) 



(0,1) 



2ms •/ 2m3 
Spin-dependent factors may be obtained with the trace formalism. One has 



(Q\qiY Ji v \h v J) . 



(fx (0\TqiT^h v 2 {x)\h v J) 



■|(-0oo 



A 2 ) tr \iT^Mj\ 



V2 



-IUJ„ 



(8.8) 



J d 4 x (OlTqiT^hv B (x)\h v J) = -i 



3 A B )tr [ir^^Mja^]^ 



dj<pooA B 



3V2 



-iuj,. 



(8.9) 



(Olgir^, h v \h v J) = - 



A 3 tr \ir^ZMjj 



d.jcb^A 



V2 



-IUJ, 



where A 2 , A B , and A 3 parametrize the light degrees of freedom, and d H 
Combining the equations of motion and heavy-quark symmetry, 



3, d 



H* 



A 3 = \(A 



m 



qji 



where m q is the mass of the light quark |2"2" |. Combining Eqs. 

A 2 dj A B ' 



2m 2 3 2m B 



(o,i)dj A 



2m H 



(8.10) 
= -1. 

(8.11) 
(8.12) 



As expected on the general grounds outlined in Sec. |V|, the rest mass does not appear. Pre- 
viously this had been shown only by explicit calculation [ 22f . Like the mass formula fl6.12|) , 
this result is simple enough that it could have been written down upon inspection of the 



corresponding continuum formula ||40| , |22 



To obtain the correct static limit of the decay constants, one must adjust the normaliza- 
tion factors Zj to yield rjj in the leading terms. This is known at the one-loop level for the 
Wilson JO] and Sheikholeslami-Wohlert actions [Q. Similarly, to obtain the 1/mg correc- 
tions, one must adjust the lattice action and currents so that m 2 = m B = = uiq, which 
is easy at the tree level. With these choices, Eq. ( |8.12| ) predicts that the heavy-light decay 
constants should depend mildly on the lattice spacing. Explicit calculation supports this 
prediction |43",44j|. On the other hand, when not all these choices are made, the dependence 
on the lattice spacing could be more pronounced, because then l/m 3 or l/m B could vary 
rapidly with m^a. Explicit calculation supports this prediction too [[4^ . 
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IX. DISCUSSION AND CONCLUSIONS 



Two themes run through Symanzik's application of effective field theory to the study of 
cutoff effects. The first is descriptive [0. The local effective Lagrangian organizes devia- 
tions from the continuum limit through a series of higher- dimension operators, multiplied 
with certain coefficients. When the higher-dimension terms are small, they can be treated 
as perturbations, and their influence can be propagated from the effective Lagrangian to 



physical quantities. The second theme turns the description into a weapon 16 . Details 
of the underlying lattice action alter the effective Lagrangian only via the short- distance 
coefficients. If a given action leads to a reduced (or vanishing) coefficient, then the process 
independence of the coefficient guarantees that its associated operator has a reduced (or 
vanishing) effect on all observables. 

The two themes also run through the application of HQET to lattice QCD. The con- 
crete results — the expansions given in Eq. ( p.1'2 ), Eqs. ( |7.2(J| )-( |7.28| ), Eqs. ( [7.3Up -( [7.34p , 
Eqs. ( |7.46| )-( [7.49| ), and Eq. flS.I2j ) — describe the deviations from the static limit of the 
mass, semileptonic form factors, and decay constant of heavy-light mesons. These descrip- 
tions hold, as always in HQET, when momentum transfers are much smaller than the heavy 
quark mass(es). Details of the lattice alter the validity of the description superficially: they 
merely change the short-distance coefficients. On the other hand, the details alter the utility 
of the description greatly: if a coefficient is tuned correctly, to some accuracy, in one observ- 
able, then its associated operator contributes correctly, to that accuracy, in all observables. 
In all examples, one sees that the leading 1/mg dependence is reproduced correctly if the 
short- distance coefficients l/rri2, 1/rriB, and l/m.3 are adjusted correctly. These conditions 
can be obtained, respectively, through suitable adjustments of the bare mass, of the "clover" 
coupling in the Sheikholeslami-Wohlert action, and of a tunable parameter in the current. 

It may be worthwhile to contrast the formalism developed here with other methods for 
treating heavy quarks in lattice gauge theory. One approach is to derive HQET or NRQCD 
in the continuum and discretize the result. In fact, both effective theories were originally 
formulated with this idea in mind |§-12]. The resulting lattice theory has ultraviolet diver- 
gences that are more severe than those of QCD, so one must either keep a -1 ~ txiq and 
employ a highly improved lattice action |fO|,|nj or restrict one's attention to the leading term 
of the infinite-mass limit p5 |. The approach developed here and in Ref. |TJ examines the 
large-mass limit of Wilson fermions, and as a — > the only ultraviolet divergences that are 
encountered are those of QCD. 

Another approach is based on lattice actions that are asymmetric under interchange of 
the temporal and spatial axes With a suitable adjustment of the asymmetric couplings, 
the physics can be made relativistically covariant. For example, one can adjust the action so 
that rrii = m>2- Cutoff effects can be analyzed either with Symanzik's effective Lagrangian, 
provided one retains the full dependence on rriQa in the coefficient functions, or with the 
HQET description developed here. Initial results [0 with the asymmetric action indicate 
that the Symanzik and HQET interpretations give the same physical results. 

Many papers have followed an ad hoc combination of Symanzik and heavy-quark effec- 
tive theories. Numerical data are generated with artificially small heavy-quark masses, to 
reduce rriQa. Then these data are extrapolated up in mass guided by the (continuum) 1/rriQ 
expansion. In practice, however, it is hard to find a region with mqa <C 1, for Symanzik's 
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analysis genuinely to apply to cutoff effects, and Aqcd/^q "C 1, for HQET genuinely to 
apply to the mass dependence. Often neither asymptotic condition realistically describes 
the numerical data. The description developed in this paper naturally applies to the subset 
of such data where HQET is indeed valid, so these data could be reanalyzed in light of the 
expansions given above. 

One might also imagine reducing the lattice spacing a by an order of magnitude or so. 
In this regime, the pictures painted by HQET and Symanzik's effective Lagrangian become 
indistinguishable from each other, even for the bottom quark The brute-force approach 
is costly, however. Processor requirements grow as a -5 (if not faster) and memory as a -4 . 
For B physics it makes more sense to invest steady improvements in computers into removing 
the quenched approximation, rather than into a radical reduction of the lattice spacing. 

Some readers may consider the proliferation of short-distance coefficients for higher- 
dimension operators to be impractical. The proliferation is genuine: it appears to the same 
extent in lattice NRQCD and to almost the same extent in the usual HQET. In many cases, 
however, the uncertainty from one-loop coefficients is smaller than other numerical uncer- 
tainties p5|j36[1 . Moreover, with limited computer resources the approach of this paper is 
more practical than a brute-force reduction of the lattice spacing and on sounder theoretical 
footing than ad hoc fitting procedures. 

A gap left by this paper is the calculation of the short-distance coefficients, which depend 
on the lattice action. Detailed calculations have been and will be addressed elsewhere |21~1 , [38| . 
The coefficients can be obtained with some accuracy via perturbation theory in the gauge 
coupling. There are, for example, general formulae, valid to every order in perturbation 
theory, relating the self energy of the underlying lattice theory to the first two coefficients of 



the effective Lagrangian, m\ and I/7B2 [21]. Similarly, radiative corrections to the currents 
are related to the (on-shell) vertex function p8|| . Beyond the one-loop level the calculations 
will not be easy, but at least they are well defined. 

An even better strategy would be to devise non-perturbative methods for tuning, if not 
explicitly calculating, the short- distance physics. For example, heavy-quark expansions of a 
hadron's kinetic mass, chromo magnetic mass, etc., would be useful, because with them one 
could remove HQET scheme dependence. Other possibilities might mimic strategies invented 
for light quarks, such as imposing — at finite lattice spacing — identities of the continuum 
limit. For heavy quarks, reparametrization invariance ])TT| , which is closely related to Lorentz 
invariance and heavy-quark symmetry, may be helpful. 

The heavy-quark expansions in this paper are just the beginning. A wide variety of phys- 
ically interesting observables have been studied with the usual HQET, and matrix elements 
of the infinite-mass limit are almost always needed. One can re-analyze each observable 
with the modified coefficients appropriate to the HQET description of lattice gauge theory, 
to find out how a direct lattice calculation compares to the continuum. Furthermore, it 
might be possible to extract parameters such as A and Ai by calculating the short- distance 
coefficients (in a suitable scheme) and fitting lattice data. The idea is similar to a pro- 
posal j|SJ for extracting kaon matrix elements from current-current correlation functions 
(J(x)J(O)). (A significant difference is that here the ratio mga of short distances is treated 
exactly, whereas in Ref. [[48] the analogous ratio a/x is presumed small.) Determinations 
of A and Ai are intriguing, because they also appear in heavy-quark expansions of inclusive 
processes [pirafM . 
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APPENDIX A: TRACES FOR SEMI-LEPTONIC FORM FACTORS 

This appendix gives the traces needed to express the semi-leptonic form factors, at zero 
recoil. Matrix elements with v' = v are considered first, in Appendix AT. They enter 



into h+(l), hi(l), and /^(l). To extract h~(l) one must take v' different from v, focus on 
terms multiplying |(V — v ) M , and then set w = 1; cf. Appendix A 2 . 



1. At zero recoil 

The traces needed to express matrix elements used to obtain h + (l), hi(l), and /^(l) 
are worked out here. One finds no contribution of the types (j^) and (j^C^) when w — 1. 

a. Contributions from (j^) 

At leading order in the heavy-quark expansion, all matrix elements are written 

{c v iJ'\c v tTb v \b v J) = -tr{Mj>rMj}£(w) (Al) 

where w = —v' ■ v. The spin dependence factors completely; there is only one function £(w) 
to parametrize the light degrees of freedom. At zero recoil the current iv^c v h v is the Noether 
current of heavy-quark flavor symmetry. The associated charge changes nothing but the 
heavy-quark flavor, namely 

d 3 y(c v J\iv°c v b v (y) = (b v J\, (A2) 

and hence £(1) = 1. Fortunately, this conclusion does not depend on the conservation of 
the current in the underlying theory, because for lattice QCD one usually computes the 
transition with a current that is not conserved. (That is why Zy is written explicitly.) The 
violation of current conservation is a short- distance effect, however, so it can appear only in 
the short- distance coefficients. 

The matrix elements of interest are 

(cvO|cv*7^M> = \{v' + vT£(w), (A3) 

(c v ,l\c v .i'fb v \b v l) = \iy' + vY7-e^(w), (A4) 

(c v ,l\c v >i^ l5 b v \b v 0} = |[(1 + w)i^ + z? • vv'^}i{w). (A5) 

(c v/ 0\c v ,i^ 5 b v \b v l) = -|[(1 +w)ie fl + ie-v'v^(w). (A6) 
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In Eqs. (|A5|) and ( |A6| ) note that e' • v = and e • v' = at zero recoil. 

The Isgur-Wise function £(u>) is ubiquitous, reappearing, for example, in h-(w), which 
is considered in the next section. Here we are concerned with v' = v , and then 

(c v J'\c v iT^b v \b v J) = = uj", (A7) 

where = 7^ or 7^75 and uj^ = i> M , v^e'-e, ^'e /M , or — ze^, as the case may be. 

b. Contributions from (j^'jC^ 1 ') 

The dimension-five interactions in the HQET Lagrangian lead to time-ordered products 
of with O2 and Ob- For unequal velocities the matrix elements are parametrized by 
three functions 

/ d A x (c vl J'\c vl TbJS)O b Jx)\b v jy = -tT{Mj>rMj}Ai(w), (A8) 

J-T 

[ d 4 x (c v ,J'\c v ,Tb v (0)O b B {x)\b v jy = -ti{Mj'Ts a(S MjA a/3 (v, v')}, (A9) 

J-T 

where, like the chromomagnetic field B a/ 3, the tensor A a p(v, v') is anti-symmetric and 
v a A a p(v,v f ) = 0. A general decomposition satisfying these constraints is 

A a p{v, v 1 ) = (r)iar)) a pA 3 (w) + (ijtav'xp ~ iv' La ^ ± p)A 2 {w). (A10) 

The same functions appear for insertions of 0\ and O c B . 

One can work out the traces to see how A\{w) and A 3 (w) contribute to h + (w), hi(w), and 
hA 1 (w). [A 2 (w) contributes to h-(w).] We are, however, mainly interested in the zero- recoil 
point, w = 1. Then the currents become Noether currents, and there are further constraints. 
With one insertion the starred time-ordered product is identical to the connected one: 

(c„J'|^r6„(0)<I^(:r)|M)* = (c v J'\c v Tb v (0)O b x (x)\b v J) c = 

(c v J'\c v Tb v (0)O b x (x)\b v J) - (^J'lcAmbv^iv^-'MO^lKJ), (All) 

for x° < 0, as in Eq. ( |A8| ) . By translation invariance the left-hand side of Eq. (|A8|) 

Jd 4 x(c v J'\Tc v Tb v (y)O b x (x)\b v jy = 1 Jdx°d 3 y(c v J'\Tc v Tb v (y)O b x (x)\b v jy. (A12) 

Taking T = iv° and using Eq. ( |A2j ) one sees that the right-hand side of Eq. ( |A11| ) vanishes 
identically. Thus, 

J d A y (c v J'\Tc v Tb v (0)O 2 (y)\b v jy = 0, (A13) 

J d A y(c v J'\Tc v Tb v (0)O B (y)\b v jy = 0, (A14) 

namely Ai(l) = and A 3 (l) = 0. 

These results are properties of heavy-quark symmetry and not of the underlying theory. 
Usually it is argued that Ai(l) = A 3 (l) = as a consequence of current conservation 
in QCD. This line of argument would not have been enough for our purposes, because for 
most choices of V^t current conservation fails. Fortunately, the foregoing argument does 
not rely on the underlying theory; indeed, it is equivalent to the derivation in Rayleigh- 
Schrodinger perturbation theory of the Ademollo-Gatto theorem. 
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c. Contributions from (j^°'C^ 2 ') 



By the same argument leading to Eqs. ( |A13| ) and (A14) 

d 4 y (c v J'\Tc v Tb v (0)O D (y)\b v jy = 0, 
A (c v J'\Tc v Tb v (0)O E (y)\b v jy = 0. 



(A15) 
(A16) 



The same holds for insertions of the four-quark operators omitted from Eq. (|3.17 ). Again, 
this is a property of heavy-quark symmetry and not of the underlying theory. 
References [0,0] choose a basis with the operator 

Qd = 2h v D^{-iv ■ V)D^h v , (A17) 

and a similar, spin-dependent operator Q E , instead of Od and Oe- They are related by 

Q d = O d + h v D 2 ± (-iv ■ V)h v + h v (iv ■ V)D 2 Ji v , (A18) 

up to total derivatives, and similarly for Qe- The additional terms, which superficially 
vanish by the equations of motion, generate contact terms. Thus, 

Jd 4 y(c v J'\Tc v Tb v (0)Q b D (y)\b v jy = (c v J'\c v TD 2 ± b v \b v J) = -tr{Mj>TMj}\i, (A19) 
Jd 4 y(c v J'\Tc v Tb v (0)Q b E (y)\b v jy = (c v J'\c v TBb v \b v J) = -ti{Mj / Ts^Mjtcr a(3 }X 2 , (A20) 

where B = s a ^B a p, and Ai and A2 are the same constants (including any light-sector cut- 
off effects) as in Eqs. ( |6.5| ) and (fT6|). The left-hand sides of Eqs. ( |A19| ) and ( |A20| ) were 
parametrized, respectively, with 2.B 1 (1) and 2B 3 (1) in Ref. ||3!| and with 2S X and 2S 3 in 
Ref. [0| , but the identification with Ai and A2 was not made. 

In the basis employing Oe and Oe, the counterpart of these contact terms are the 
contributions c v >TDj_b v and c v >TBb v to the currents, cf. Eqs. ( [7.16[ ) and ( [7.17[ ). In the Q- 



basis these currents have coefficients (8m^ 2 ) 1 — (8m|,) 1 and {%m 2 sB ) 1 — (8m E ) 
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d. Contributions from (j^) 

There are two kinds of of second-order corrections: those which can be associated with 
a single leg and those which involve cross-talk between the legs. At zero recoil all can be 
expressed through the parameters Ai and A2, namely 

(c v J'\c v rV a V\\b v J) = -tr [Mj'TMj [iAir? Q/3 + \\ 2 ia aP ] } . (A21) 

By taking V to be the unit matrix or s a p and contracting indices, it is easy to trace back to 
the definitions ( |6.5| ) and (|6.6|) . By dimensional analysis, these are the only corrections that 
can arise, even beyond tree level. 
The required matrix elements are 
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{c v J'\c v iT ^b^bvJ) = {c v J'\c v D JT ^b v \b v J) = Ai^, (A22) 

(c v J'\c v iT^Bb v \b v J) = dj\ 2 u}^ (A23) 
(c v J'\c v BiT^b v \b v J) = dji\ 2 uj^ (A24) 

At zero recoil (c v J'\c v V a TV%\b v J) = -(c v J'\c v TV a V%\b v J), so and have 

matrix elements 

(cvJ^p^^p L b v \b v J) = (Ai + dj\ 2 )u}^ (A25) 

(c u l|c u |>±i7At75-Pj.6i>|6«0) = (c„0|c, ) |)_ L z7 A1 7 5 |)_ L 6„|6, ) l) 

= -|(Ai + 3A 2 )^ (A26) 

Contributions with Ai (A2) are spin-independent (spin-dependent). 

e. Contributions from {jO^j^C^) 

Several matrix elements are introduced for double insertions of In the following 

the short- distance coefficients are stripped off, leading to insertions of / d 4 zO\(z) 1 where 
X G {2, B} and h labels the heavy flavor. When the operator comes from the numerator 
of Eq. (|5.4j ) the time variable is integrated for h = b over the interval (— T, 0] and for h = c 
over [0, T); when the operator comes from the denominator the time variable is integrated 
over the interval (-T, T). After generating all terms the limit T — > oo(l — i0 + ) is taken. 

When two interactions occur on the incoming line 

\ J d 4 xd 4 y (c v J'\Tc v tTM0)O b 2 (x)O b 2 (y)\b v jy = u^A (A27) 

J d 4 xd 4 y (c v J'\Tc v iTMQ)O b 2(x)O b B (y)\bvJy = uj^djB (A28) 

\ J d 4 xd 4 y (cJ'\Tc v tt^)0 B {x)O h B {y)\b v jy = ^(Ci + djC 3 ) (A29) 
and, similarly, when two interactions occur on the outgoing line 

I / d 4 xd 4 y {c v J'\TO c 2 {x)O c 2 {y)c v iTM^JY = uj^A (A30) 

J d 4 xd 4 y (c v J'\TO c 2 (x)O c B (y)c v iTMQ)\bvJy = ^dj,B (A31) 

\ J d 4 xd 4 y (c v J'\TO c B (x)O c B (y)c v tT ll b v mbvJy = w»{C x + d r C z ) (A32) 
where T M = 7^ or 7^75, as the case may be. When each line has one interaction 

J d 4 xd 4 y (c v J'\TO c 2 (x)c v iTMQ)£>2(y)\bvJy = oj^D (A33) 
J d 4 xd 4 y (c v J'\TO c 2 (x)c v iTM )°B(y)\bvJy = u^djE (A34) 
J d 4 xd 4 y (cJ^TO B {x)c v iT^MO b 2 {y)\b v JY = u^djrE (A35) 
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again where T M = 7 M or 7^75, as the case may be, and 

J d i xd i y(c v J\TO%(x)c v i ltM b v ( )°B(y)\bvJ)* = uj^Ri + djRi), (A36) 
J d A xd A y (c v l\TO%(x)c v i^ l5 b(0)O b B (y)\b v 0y = -\u)p(Ri + 3i? 2 ). (A37) 

Here we have used the notation of Ref. [[34]]. 

There are relations between these parameters, which follow solely from heavy-quark 
symmetry and properties of perturbation theory. Upon expanding Eq. (|5.4j ) and sorting 
terms with like coefficients one finds 

j d A xd A y(c v J'\Tj(0)O b x (x)O b Y (y)\b v jy = 

J d A xd A y (c v J'\Tj(0)O b x (x)O b Y (y)\b v J) c - (c v J'\j{0)\b v J)Z XY , (A38) 

and 

J d A xd A y (c v J'\O x (x)j(0)O b Y (y)\b v jy = J d A xd A y (c v J'\O x (x)j(0)O b Y (y)\b v J) c , (A39) 

with limits of integration on the time coordinates as given above. On the right-hand side of 
Eq. ( |A38| ) the second term is the contribution from state renormalization: 

Zxy = 4 [ T d 4 x f° d A y(b v J\O b x (x)O b Y (y)\b v J) c , (A40) 

V u JO J-T 

which is flavor independent. In Eq. ( |A38| ) the operator j is left-most for all time orderings. 
When j is a Noether charge one can apply Eq. ( |A12| ) to show that the connected term 
vanishes, leaving only the term from state renormalization. In Eq. ( |A39|) the operator j is 
in the middle for all time orderings. When j is a Noether charge, however, the right-hand 
side can be reduced to the same quantity as in the state renormalization. Inserting complete 
sets of states on both sides of j, and noting that Eq. ( |A12| ) applies equally well to excited 
states, one finds 

J d A xd A y (c v J'\O x (x)j(0)O b Y (y)\b v jy = (c v J'\j(0)\b v J)Z XY , (A41) 

making use of the flavor independence of Z XY - Apart from a sign, therefore, the two kinds 
of ^-products are the same, and 

A = -ID, (A42) 
B = —E, (A43) 
Cx = -\Ri, (A44) 
C 3 = -\R 2 . (A45) 

These identities leave only four parameters. To my knowledge they have not been derived be- 
fore. Since they do not depend on the underlying theory, they hold also for continuum QCD. 
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2. Near zero recoil: h-(l) 



At zero recoil several matrix elements vanish, but they are precisely of the type leading 
to h__(w) in Eq. (|7.2|) . To extract h_(l) one must take v' ^ v while evaluating matrix 
elements, read off the form factor, and then set w to 1. This subsection works out the 
relevant matrix elements, those of the dimension-four currents, the dimension-five currents 
c v ii^^i$b v and c v 4p i^^K, and time-ordered products of dimension-four currents with 



a. Contributions from (j^) 

For the matrix elements (c v >0\c v >i'y l j j p ± b v \b v Q) and (c v iO\c v 'p ^ij^bvlbvO) one starts with 
the matrix element 

(c v ,J'\c v ,TV a b v \b v J) = -Xx{Mj>VMM a {v,v')} (A46) 

where £ Q parametrizes the light degrees of freedom. The equation of motion (—iv ■ T>)b v = 
implies that v a £ a (v,v') = 0, leaving two independent form factors 

?(v, v') = v'?Z 2 (w) - i^lUw). (A47) 
A further constraint on £ a (i>, v') comes from the "integration- by-parts" identity 

(c vl J'\c v tV a Tb v \b v J) + (c v ,J'\c v ,TV a b v \b v J) = -iA(v' - v) a (c v r J'\c v ,Tb v \b v J), (A48) 
where T>b v — (D — imif,v)b v and c v <T> = c v i(D + imi c v'). The first matrix element 



(c v ,J'\c v ,V a rb v \b v J) = - tT{Mj>rMj[-i£ a {v', v)}}, (A49) 



where £"(V, v) = 7 4 [£ a (t> w)]^ 4 . Substituting traces for matrix elements in Eq. ( |A48| ) yields 
the relation 

(w + lMw) + &(w) = -MM, (A50) 



which can be used to eliminate ^{w). In Eq. ( |A50|) the constant A and the function £( 



[w 



are the same — including lattice artifacts of the light degrees of freedom — as in Eqs. ( |6.3| 
and ( |A1| ), respectively. Evaluating the traces of interest and using Eq. ( |A50| ) one finds 

(c v >0\c v 4'y fl p J _b v \b v 0) = (cv0| M) 



i f„.' 



v'-v)^ s (w)-M(w)], (A51) 

There is no contribution to h + (w), and the vector-to- vector matrix elements make no con- 
tribution to hi(w), just to other form factors that are not considered in this paper. An 
equivalent analysis appears in Ref. ]33|]. The only significant addition is to extend to lattice 
QCD the identification of A£(iu) in Eq. ( |A50|) with the quantities in Eqs. (|oT3| ) and (|A1|) . 
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b. Contributions from (j*- 2 **) 

To obtain all of the second-order corrections to the current one can start with 

(c v ,J'\c v ,V a TVPb v \b v J) = -tT{M r TMj[-\ aP (v,v')]}. (A52) 

The equations of motion {—iv ■ T>)b v = and c v >(iv' ■ T>) = imply that A Q/3 (t> , v ')vp = 
and v' a X al3 (v , v') = 0, and symmetry under exchanging final and initial states implies that 
\P a (v',v) = \ a/3 (v,v'), leaving four independent form factors, 

\<#(v, v') = V '^ s X 1 (w) + \ia^\ 2 {w)}4 + v a ± ,v'{\ 3 (w) + [i^v' P ± + v^i-y^w). (A53) 

The pre-factors for the first two form factors are chosen so that Ai(l) = Ai and A2(l) = A2 
are the constants in Eq. ( |A2 1| ) . 



The matrix elements needed for h-(w) are (c v i0\c v >i'y l j i i$b v \b v 0) and (c v '0\c v >ip'i'j fl b v \b v 0). 
They are related to Eq. ( |A52| ) by the identity 

(c v ,J'\c v ,TV a V\\b v J) = -iA(v' - v) a (c v ,J'\c v ,TV p b v \b v J) - (c v ,J'\c v ,V a TV\\b v J} 

(A54) 

and the definitions if* = —iv a [V a ,Vf 3 ]'y±, if!' = —iv^lV ,!)^}^,. Evaluating the traces one 
finds 

(c v '0\c v 'i , y IJ 'if]b v \b v 0) = (c v >0\c v >ip'i'y >1 b v \b v 0) 

= -\{v' - vY{\{w + l)A(tu) + {w- 1)A [2£ 3 {w) - M(w)] } (A55) 

where 

\(w) = IwX^w) + (3 - w)\ 2 {w) - 2(w 2 - 1)A 3 (V) + 8(w - l)A 4 (u;). (A56) 
Atw = l, A(l) = |(Ai + 3A 2 ). 

c. Contributions from (j^C^) 

The time-ordered products of interest are / d 4 y (c v > J'\T c v /Tp ± b v (x)Ox(y)\b v J)* and 

fd 4 y(c v ,J'\Tc v/ p ±/ rb v {x)O f x {y)\b v J)\ where X e {2,B} and / G {c,b}. As before it 
is helpful to consider matrix elements with p ± replaced with T> and derive constraints from 
the equations of motion and from "integrating by parts." This is a bit trickier now, with 
derivatives acting under the time-ordered product. 
The equations of motion imply the identities 

J d A y (c v ,J'\Tc v/ T(-iv ■ V)b v (x)O b 2 (y)\b v jy = (c v > J' | c v >YD 2 A b v (x) \b v J), (A57) 

J d A y (c v/ J'\Tc vl T(-iv ■ V)b v (x)O b B (y)\b v jy = (c v > J'\c v ,TBb v (x)\b v J) , (A58) 

d 4 y(c v/ J'\TO c 2 (y)c v ,(iv' ■ V)Tb v (x)\b v J)* = (c v >J'\c v >D 2 1 _ l rb v (x)\b v J), (A59) 

d 4 y (c v/ J'\TO c B (y)c v ,(iv' ■ V)Tb v (x)\b v jy = (c v/ J'\c v >B'Tb v (x)\b v J}. (A60) 
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The contact terms on the right-hand side were omitted from Eqs. (4.27) of Ref. [3^] but do 



appear, for example, in Eq. (A21) of Ref. [49|. They arise from a careful definition of the 
T-product for operators containing time derivatives. A helpful mnemonic for checking them 
is to note that 



{-ivV)Tb v {x)b v {y)=5^\x-y) y 
Tc v/ (y)c v ,(x)(iv' -V) =5 {A \y-x). 



(A61) 
(A62) 



Further identities come from taking the derivative d p (D v >\c v 'Tb v \B v ) between fully dressed 
states, and generating the expansion. This leads to 



d A y (c v ,J'\T [c v ,TV p b v (x) + c v ,V p Tb v (x)]O b 2 (y) 



b v J)* — i\\v p (c v i J 1 \c v iTb v (x)\b v J) , 
b v J}* = (A64) 
b v J)* - idj\ 2 v p (c v iJ'\c v <Tb v (x)\b v J}, 
b v J)* = (A65) 
b v jy + i\ l v' p (c v i J'\c v iTb v (x)\b v J) , 
b v jy = (A66) 
b v J)* + idji\2V lp (c v 'J'\c v 'Tb v (x)\b v J). 

These identities do not agree with analogous ones from combining Eqs. (C4) and (C5) 



-iA{v' -v) p I d 4 y (c v ,J'\T c v ,Tb v (x)O b 2 (y) 
d 4 y {c v ,J'\T [c v ,TV p b v (x) + c v ,V p Tb v (x)]O b B (y) 

-iA(v' - v) p J d*y (c v/ J'\Tc v/ Tb v (x)O b B (y) 
d A y (c v ,J'\T O c 2 Xy)[c v ,TV p b v (x) + c v 3 p Tb v (x)} 

-iA(v'-v) p J d A y (<vJ'|T O c 2 (y)c v ,Yb v (x) 
d 4 y {c v ,J'\TO B (y)[c v ,TV p b v (x) + c v 3 p Tb v (x)\ 

-iA(v'-v) p [ d A y(c v ,J'\TO B {y)c v ,Tb v {x) 



b v JY 



(A63) 



of Ref. |33[]. Remarkably, Eq. (C5) of Ref. [33] contains the contact terms omitted from 
Eq. (4.27) of Ref. |§. 

Once again the time-ordered products are parametrized by form factors. Consider first 
the case with the kinetic operator. It is enough to present the details for O b 2 . One may write 

J d*y (c v/ J'\Tc v/ TV a b v O b (y)\b v jy = -tr{Mj,TMjiZ a (v,v')}, (A67) 
d 4 y (c vl J'\Tc v S a rb v O b (y)\b v jy = -^{M JI YMj[-iF a {v,v')\}. (A68) 

In the first case, the equation of motion ( |A57| ) implies v ■ H = <po, where 

0o H = |(2 + w 2 )AiH -(w- l){w[±\ 2 (w) + (w + l)X 3 (w) - 2X 4 (w)} + A 2 £(w)} (A69) 
is obtained from Eqs. ( |A53| ) and (|A54J ) . Thus, S a has a decomposition 

' -i;>oH+<S 2 H-z 7 ^3H (A70) 



iv,v 



similar to £ a but with —<f) (w) multiplying v a . On the other hand, the equation of motion 
still implies v' ■ F — 0, so F a has the decomposition 
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F a (v, v') = vl,F 1 (w) - i"ft,F 3 (w). (A71) 



similar to £ a (v',v). The form factors 5 2 , Fi, and F 3 can be eliminated, because the iden- 
tity ( |A63| ) implies 

(w + 1)S 2 + ~ 3 = -w0o - AAi, (A72) 
(w + l)Fi + F 3 = 0o - AAi, (A73) 
^3 = S 3 , (A74) 

where 

Mw) = frM- A ^) . (A75) 
ty — 1 

At zero recoil the new constants that can arise are S 3 (l) and, denoting differentiation with 
respect to w by a dot, </>o(l) = 0o(l) — Ai^(l). (Recall that Ai(l) = 0, as a consequence of 
heavy-quark flavor symmetry.) 

Evaluating the traces for h-(w), one finds 



d 4 y (c v ,0\Tc v n^p x b v (x)O»(y)\b v 0y = J d 4 y (c v ,0\Tc v ,O c 2 (y)p ± aYb v (x)\b v 0)* 

= -§(t/ - vY[2E 3 (w) - AA 1 (w) - wMw)}, (A76) 
d A y (c v/ 0\Tc v ,p ± nYb v (x)O b 2 (y)\b v 0y = J d A y (c v ,0\T c v ,O 2 (y)i-f p p ± b v (x)\b v 0y 

= -\{v' - <[2~ 3 H - AA 1 (w) + Mw)}. (A77) 

Matrix elements of this kind make no contribution to h+(w), hi(w), or Iia^w). 

Finally there are the time-ordered products with the chromomagnetic energy. It is enough 
to show the details for O b B . When the derivative acts on b v , 

J d 4 y (c vl J'\Tc v ,TV p b v {x)O h B {y)\b v jy = -tx{Mj,Ys^ ''MjiZ p a p(v,v')}, (A78) 

and when the derivative acts on c„/, 

J d 4 y (c v ,J'\Tc v ,V<>Tb v {x)0 B {y)\b v jy = -tx{Mj>Ys^Mj[-iF' , a p{v,v')\}. (A79) 

The tensors 'Ef a p and F p a p inherit properties from the chromomagnetic field B a p\ they are 
antisymmetric on the lower indices, and S p Q/ gf /3 = F p a pv" = 0. From the equations of 
motion v' p F p a p = and v p E p al 3 = <f> a p, where 

Mv, v') = ^[^{v, v') - \"»(v, v')]r} v p + AKJp(v, v') - v' ±p Uv, v% (A80) 
Substituting Eq. ( |A5g ) into Eq. (TA"80| ) 

<f> a p{v, v') = {r}icrr]) a p(j) 3 {w) - (ij± a v' ±j3 - iv' ±c jx/3)(f>2(w), (A81) 

where 3 (w) = \2(w) and <f>2(w) = — h\ 2 (w) — (w + l)A 4 (w) — A£ 3 (w). The constraints 
on E p a/ 3 and F p a p lead to the decompositions 
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2 p a/3(>, v') = (r)iar)) a/3 [-v p (j) 3 + v'[E s - i^Eg 

- (il-L a v'±p - v' La ii±p) [-v p <p 2 + v'[E 5 + r^H 6 ] (A82) 
+ (Va v ±/3 - Vpv' ±a )E 10 + (r)&j ±/3 - ^i7± a )Sn, 



and 



F p af 3(v, v') = {r]icrr}) a p [v p ± ,F 7 - Z7^,F 9 ] 

- (hxaVxp - v' ±a i^±p) [v p ±l F 4 + i^j_,F 6 ] {A83) 

The subscripts are chosen as in Ref. [Q. 

The identity ( |A64| ) can be applied to eliminate S 5 , E 8 , and all Fs: 

F k = E k , k G {6, 9, 10, 11}, (A84) 

(w 2 - 1)~ 5 = -wfo ~(w + 1)S 6 - 5n + (w- 1)AA 2 , (A85) 

(w 2 - l)F A = 2 + (w + l)Se + wE n + {w- l)AA 2 , (A86) 

(w + 1)S 8 = -u>0 3 - S 9 - AA 3 , (A87) 

(w + 1)F 7 = 03 - E 9 - Av4 3 , (A88) 

where 

= ^sH-W (Agg) 

10 — 1 

Each of Eqs. (fAS§) and flABBp implies 2E 6 (1) + H u (l) = -0 2 (1). 
Evaluating the traces for h-(w), one finds 

d 4 y (c v ,0\T c v ,ij p p ± b v (x)O b B (y)\b v 0y = J d 4 y {c v ^T c vl O c B {y)p ^ p b v {x)\b^Y 

= -\{v' -vY[2E_(w) - w4>-(w)], (A90) 
d A y (c v ,0\Tc vl p ± a^b v {x)O b B {y)\b v 0)* = J d 4 y (c v ,Q\T c v ,O B {y)i^p x b v {x)\b v 0)* 

= -|(v / -v)'*[2S_H + 0_(w)], (A91) 

where 

E = 3E 9 + (w + 1)(2S 6 + S 10 ) - 2S U - §AA 3 - (w - 1)AA 2 . (A92) 
0_ = 303 - 20 2 (A93) 

At zero recoil the new constants that can arise are ^(l), S 9 (l), Sio(l), and 03(1). (Note 
that A 2 (l) drops out, and recall that ^3(1) = as a consequence of heavy-quark flavor 
symmetry.) As in Eqs. (|A76|) and ( |A77| ), matrix elements of this kind make no contribution 
to h+(w), h\(w), or Iia^w). In h-(l) they reduce to two constants 

S_(l) = 359(1) + 8S 6 (1) + 2Sio(l) + 20 2 (1), (A94) 
0_(1) = 3[0 3 (1) - A 2 £(l)] - 20 2 (1), (A95) 

which are needed in Eq. ( [7.47] ) . 
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